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SUMMARY 


The present report presents a theoretical analysis of an initially 
flat, rectangular plate with largo deflections under either normal 
pressure or combined normal pressure and side thrust. As small 
deflections of a flat plate are governed by a single linear equation, 
large deflections introduce nonlinear terms into the conditions of 
equilibrium and are governed by 2 fourth -order, second -degree, partial 
differential equations. These so-called Von Karman equations are 
studied in the present report by use of the finite -difference approxi- 
mations. The difference equations are solved by two methods, namely, 
the method of successive approximations and the relaxation method. 
Neither of these methods is new, but their application to nonlinear 
problems requires new techniques. 

The problem of a uniformly loaded square plate with boundary 
conditions which approximate the riveted sheet-stringer panels is 
solved by the method of successive approximations. The theoretical 
center deflections show good agreement with the recent experimental 
results obtained at the California Institute of Technology when the 
deflections are of the order of the plate thickness. This agreement 
perhaps suggests the range in which these Von Karxoan equations are 
to be applied. 

Other problems of thin plates with large deflections are discussed 
from the point of view of an aeronautical engineer. The boundary 
conditions which approximate the various cases are formulated, and the 
methods for solving these problems are outlined. 

Since the method presented in the present report is general, it 
may be applied to solve bending and combined bending and buckling 
problems with practically any boundary conditions, and the results may 
be obtained to any degree of accuracy required. Furthermore, the same 
method may be applied to solve the membrane theory of the plate which 
applies when the deflection is very large in comparison with the 
thickness of the plate. 
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INTRODUCTION 


The classical theory of the "bending of a thin elastic plate expresses 
the relation between the transverse deflection of the middle surface of 
the plate w and the lateral loading of intensity p by the equation 


D V^v 


= P 


where D = — — is the flexural rigidity of the plate . It is 

12(1 - n 2 ) 

known that the theory is restricted in application, for on the one hand 
its basic assumptions can be questioned unless the plate is thin, and 
on the other hand it neglects an effect which must be appreciable 
when w has values comparable with the thickness. This is the membrane 
effect of curvature, whereby tension or compression in the middle 
surface tends to oppose or to reinforce p. The effect is negligible 
when w Is very small, provided no stresses act initially in the plane 
of the middle surface j but even so, it operates when w is small 
because stretching the middle surface is a necessary consequence of the 
transverse deflection. When the deflection gets larger and larger, 
the membrane effect becomes more and more prominent until for very large 
values of w the membrane effect is predominant whereas the bending 
stiffness is comparatively negligible. 

Small transverse displacements of a flat elastic plate are governed 
by a single linear equation but large displacements entail stretching 
of the middle surface and consequent tensions which. Interacting with 
the curvatures, introduce nonlinear terms into the conditions of 
equilibrium and so make those equations no longer independent. 

The large -deflection theory of flat plates is given by A. Foppl 
(reference l), and the second-order terms were formulated by Theodore 
von Karman in 1910 (reference 2). The amended (large -deflection) 
equations have been solved, however, in only a few cases (references 3 
to 19 ) and then with considerable labor. 

Essentially there are three problems concerning flat plates with 
large deflections. They are: 

1. The bending problems, when the flat plates are subjected to 
lateral loading perpendicular to the plane of the plates, but no 
side thrust is applied in the plane of the plates 

2. The buckling problems, when the plates are subjected to side 
thrusts In the plane of the plates but are not loaded laterally 

3. The combined bending and buckling problems, when the plates 
are subjected to both lateral loading and 3ide thrusts 
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In the case of metal airplanes, in which weight is of primary 
importance, the metal sheets used must "be thin and the deflections of 
the plates are usually large in comparison with their thickness. In 
order to obtain the design formulas or charts for stressing such 
plates, the large -deflection theory must be used. 

The bending problem is important in the design of seaplanes. 
Seaplanes are subjected to a severe impact during landing and take- 
off, especi ally on rough water. The impact must be withstood first 
by the bottom plating and then by a system of transverse and longi- 
tudinal members to which the bottom plating is attached, before it is 
transmitted to the body of the structure. The bottom should be strong 
enough not to washboard permanently under these impact pressures. 

Such washboarding is undesirable because of the increased friction 
between the float bottom and the water and also because of the increased 
aerodynamic drag in flight. 

The bottom plating of seaplanes is, as a rule, subdivided into 
a large number of nearly rectangular areas by the transverse and 
longitudinal supporting ribs. Each of these areas behaves substantially 
like a rectangular plate under normal pressure. Bending of rectangular 
flat plates may therefore be used to study the washboarding of seaplane 
bottoms, provided the boundary conditions at the edges can be formulated 
Just as in the seaplane . 

The buckling problem is important in determining the strength 
of sheet-stringer panels in end compression. The use of stiffened 
sheet to carry compressive loads is increasingly popular in box beams 
for airp lan e wings and in other types of semi monocoque construction. 
Inasmuch as the sheets used as aircraft structural elements are generally 
quite thin, the buckling stresses of these sheet elements are necessari y 
low. The designer is therefore confronted with the problem of using 
sheet metal in the buckled or wave state and of determining the stress 
distribution and allowable stresses in such buckled plates. 


The combined bending and buckling problem ha3 become a problem of 
importance with the increasing use of wings of the 3 tressed-skin type 
and the pressurized fuselage construction for high-altitude flight. 

During flight the wing is subjected to a pressure difference between 
the two sides which produces the lift. The normal pressure acts 
directly on the sheet covering and is then distributed to ribs and 
spars. At the same time tho sheet panels are also subjected to a side 
thrust due to bending of the wing. In an airplane of pressurized 
fuselage construction an attempt is made to keep tho pressure inside 
the cabin at a comfortable level for the passengers, regardless of the 
altitude of the airplane. Thus, there is a pressure differential across 
the fuselage skin with an internal pressure higher than that outside. 

The fuselage skin is usually subdivided into a number of rectangular 
curved panels by longitudinal stringers and rings. These panels are 
sub lected to the pressure difference and side thrust resulting from 
tending of the f useless . As pointed out ty Niles and Newell (reference 20) 
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the strength of curved sheet -stringer panels can he determined approxi- 
mately from the flat sheet -stringer panels. The problem is then 
essentially that of determining the strength of flat plates under combined 
lateral loading and side thrust. 


Levy (reference 19) has shown that the effective width of a square 
plate with simply supported edges decreases with the addition of lateral 


pressure and that the reduction is appreciable for 


BS- 

Eh 4 


> 2.25* Therefore. 


a panel is unsafe if its design is based upon the side -thrust considerations 
only , and the study of combined loading is of great signif icance . 


A great number of authors have studied the buckling problems, and 
considerable experimental work has been carried out. As a result, 
design formulas are available and seem to be accurate for most practical 
purposes. The bending problems, however, have been studied by only a 
few investigators, and test results (references 21 to 23) are far too 
scarce to justify any conclusions. The combined bending and buckling 
problem has been studied in only one case (reference 19) , and even in 
this instance the results are incomplete. 


Among the solutions of the large -deflection problems of rectangular 
plates under bending or combined bending and compression. Levy's solutions 
are the only ones of a theoretically exact nature. His solutions are, 
however, limited to a few boundary conditions and the numerical results 
can be obtained only alter great labor. 

The purpose of the present investigation is to develop a simple 
and yet sufficiently accurate method for the solution of the bending and 
the combined bending and buckling problems for engineering purposes, 
and this is accomplished by means of the finite -difference approximations. 

Solving the partial differential equations by finite -difference 
equations has been accomplished previously. Solving the resulting 
difference equations, however, is still a problem. In the case of 
linear difference equations, solutions by successive approximation 
are always convergent and the work is only tedious. Besides, Southwell's 
relaxation method may be applied without too much trouble. Eut, in 
order to solve the nonlinear difference equations, the successive - 
approximation method cannot always be relied on because it does not 
always give a convergent solution. The relaxation method, since it is 
nothing but intelligent guessing, can be applied in only a few cases 
and then with great difficulties (reference l6). 

A study of the finite -difference expressions of the large -deflection 
theory reveals that a technique can be developed by moans of which the 
system of nonlinear difference equations can be solved with rapid 
convergence by successive approximation by using Crout's method of 
solving a system of linear simultaneous equations (reference 24) . By 
way of illustration, a square plate under uniform normal pressure will] 
boundary conditions approximating the riveted sheet-stringer panel 
is studied by this method. Nondimens ional deflections and stresses are 
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given under various normal pressures. The results are consistent with 
Levy's approximate numerical solution for ideal, simple supported plates 
(reference 19) and Way's approximate solution for ideal clamped edges 
(reference 15 ), and the center deflections check closely with the test 

results hy Head and Sechler (reference 23) for the ratio pa /Eh as 
large as 120. The deviation for the ratio pa 4 /Eh larger than 120 

is probably due to the approximations employed in the derivation of 
the basic differential equation. 

The procedure is quite general j it may be applied to solve the 
problems of rectangular plates of any length-width ratio with various 
boundary conditions under either normal pressure or combined normal 
pressure and side thrust . 

The present investigation was originally carried out under the 
direction of Professor Joseph S. Newell at the Daniel Guggenhe m 
Aeronautical Laboratory of the Massachusetts Institute of Technology 
and was completed at Brown University, under the sponsorship and with 
the financial support of the National Advisory Committee for Aeronautics, 
where the author was participating in the program for Advanced Instruction 
and Research in Mechanics. The author was particularly fortunate to 
receive frequent advice while working on this problem from Professor 
Richard von Mises of Harvard University. The author is grateful to both 
Professor Newell and Professor von Mises for their many valuable 
suggestions . 


SYMB0IB 


a, b 
h 

x, y, z 

U, V 


w 


p 


E 


D 


length and width of plate, respectively 

thickness of plate 

coordinates of a point in plate 

horizontal displacements of points in middle surface 

in x- and y-directions, respectively (nondimensiona. 

forms are ua/h 2 , va/h 2 , respectively) 

deflection of middle surface from its initial plane 
( nondimensional form is w/h) 

normal load on plate per unit area (nondimensional 
form is pa^/Eh^) 

Young's modulus and Poisson’s ratio, respectively 

r s h 3 ^ 

flexural rigidity of plate V. pg(l - \sP\J 
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y2 . S- t 

dx 2 ^y 2 

V 1 * , 4i + a» + 

dx^ dx 2 


a ' rT 1 t * 
u x > °y > T xy 


e " c " y 1 

€ x > € y > 7 xy 


F 


A, A 2 , ..., A n 


membrane stresses In middle surface (nondimensional 
forms are a z f a 2 /Eh 2 , a'a^/Eh 2 , and T *a 2 /Eh 2 , 

j Ay 


respectively) 

°x "* CT y ", T xy " extreme -f iber bending and shearing stresses 


(nondimensional forms axe a x "a^/Eh 2 , a "a 2 /Eh 2 , 


and t "a 2 /Eh 2 , respectively) 

Ay 


€ x'> c y ' > 7 xy' membrane strains in middle surface (nondimensional 

forms are € x 'a 2 /h 2 , € 'a 2 /h 2 , and 7 * a 2 /h 2 , 


respectively) 

extreme -fiber bending and shearing strains 

(nondimensional forms are € x "a^/h2, e "a 2 /h^, 
and 7 xy "a 2 /h 2 , respectively) 

stress function (nondimensional form ie F/Eh^) 

first-, second-, ..., to nth-order differences, 
respectively 




first-order differences in x- and y -directions, 
respectively 


FUNDAMENTAL DIFFERENTIAL EQUATIONS 


The thickness of the plate is assumed small compared with its other 
dimensions. The middle plane of the plate is taken to coincide with the 
xy -plane of the coordinate system and to be a plane of elastic symmetry. 
After bending, the points of the middle plane are displaced and lie 
on some surface which is called the middle surface of the plate. The 
displacement of a point of the middle plane in the direction of 
the z-axis w is called the deflection of the given point of the plate. 

Consider the case in which the deflections are large in comparison 
with the thickness of the plate but, at the same time, are small enough 
to justify the following assumptions: 
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1. Lines normal to the middle surface before defoliation remain 
normal to the middle surface after deformation. 

2. The normal stress a z perpendicular to the faces of the plate 
is negligible in comparison with the other normal stresses. 

In order to investigate the state of strain in a bent plate, it. 
is supposed that the middle surface is actually deformed and that 
the deflections are no longer small in comparison with the thickness 
of the plate but are still small as compared with the other dimensions. 

Under these assumptions, the following fundamental partial 
differential equations governing the deformation of thin plates can be 
derived from the compatibility and equilibrium conditions: 


+ 2 -& 

af* ' ° 


- + = e 

&x 2 dy c dy^ 


<3 2 w^ 2 _ c£w c£w 
dx 2 dy 2 



F 


a 2 


s 4 w .. a 4 w . A = p + + c£f <£* _ 

+ ^c)x 2 dy 2 oy4 D - 0 vox 2 by 2 oy 2 Sx 2 : ^ x 


where D = 


Eh : 


12(1 - n 2 )’ 


the median -fiber stresses are 


and the median-fiber strains are 


V - 



a ' = 

tl 


y 

CvJ 


t _ 

B 2 f 


xy 

& 

.3 are 



if 



E\ 

_ 

^y 2 

^ % n / 

, 1 f 

^F . 

o 2 f\ 

-A 

..vJX 2 
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y . B 2(1 1 u) 5 2 F 

E dx dy 


The extreme -fiber 'bending a nd shearing stresses are 


xy 


Eh 

w 

2(1 - n 2 ) 

Vdx2 

Eh 

/X 2 w 

2(1 - n 2 ) 

w 

Eh 

c^Tr 


2(1 + ji) dx dy 
These expre ssions can be made nondimensional by writing 


F 1 = -§- 
h 2 E 


i w 
w' = - 


P* 

1 

Eh 


x' = * 


*'-S 


c =el 


where a is the smaller side of the rectangular plate . 
The differential equations then becoias 

dS" { ctX 1 \ 2 b 2 v t d 2 *' 


+ 2 jil> + 


dx f cfac' 2 cty \dx' dyy cbc' 2 dy 


i2 


a 4 - 


r + 2 + r 

ctoc 1 ^ Sx l2 dy ' 2 dy 1 ^ 


^ - 12(1 - „ 2 )p' + 12(1 - f )(\ X & 


\^y dx 


8 2 F' dX ' d 2 F' d 2 w' 


cbc 8y' 2 ox 1 $y 1 cbc ' dy> 
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Tf u 2 = o.l, which Talus is characteristic of aluminum alloys, and 

the primes a^e dropped, the partial differential equations in nondimensional 


form are 


& + 2 ii— + ti 

d? dx 2 dy 2 dy 1 * 


5 2 w'\ ^ d 2 v c> 2 v 


dyy dx 2 


3 - W + 2 , %2-S * |jj[ - 10 -^ + 10 - 8 (^2 


i\ __ 

a? " ‘a* 2 Sy 2 


(i) 


d 2 F d 2 * 3% 

dx 2 dy 2 dx dy &x dy, 


( 2 ) 


The noMimeneional median-fiber stresses are 

1 3y 2 

. 

° y = dx 2 

T _^2_ 
^ dx dy 


> 


J 


(3) 


and the nondimensional median-fiber strains are 




y " 


xy 


d 2 F 

d 2 F 

&y 2 

^cix 2 

d 2 F 

d 2 F 

dx 2 

V 

-2(1 

+ ) d 2 * 1 

+ ^3x 




J 


( 4 ) 


The nondimensional ex tren» -fiber bending and shearing stresses are 
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2(1 - & 

\chc^ 

1 

l^d 2 w 

2(1 - p 2 ) 


1 

d 2 * 


+ ^ ) > 


2(1 + ji) cbc by 


FORMULATION OF BOUNDARY CONDITIONS 


The governing differential equations are 2 fourth-order simul- 
taneous partial differential equations in two variables. In order to 
o tain a unique solution in the case of rectangular plates, there must 
he four given boundary conditions at each edge . 

Before proceeding to tho actual case, two theoretical boundary 
conditions may be mentioned: d 

rotate plates ' that ls > P lates ha ^:Lng edges that can 

rotate freely about the supports and can move freely along the supports 

C1 f iP !l^ built-in plates, that is, plates having edges that 
e clamped rigidly against rotation about the supports and at the same 
time are prevented from having any displacements along the supports 

Actually it is to be expected that neither of these conditions will 
be iulxilled exactly in a structure. 

, . The pending problem will be considered next, in which the bottom 
plating of a seaplane is to bo studied. The behavior of the sheet 
approximates that of an infinite shoot supported on a homogeneous 
elastic network with rectangular fields of the same rigidity as the 
supporting framework of the seap lan e . 

Because of the symmetry of the rectangular fields, the displacement 
in the plane of the sheet and the slope of the sheet relative to the 
plane of the network must be zero wherever the sheet masses over the 
center line of each supporting beam. Each rectangular field will 
therefore behave as a rectangular plate clamped along its four edges on 
supports that are rigid enough in the plane of the sheet to erevent 
their displacement in that plane. At the same time these sunuorts must 
have a rigidity normal to the plane of the sheet equal to that of 
the actual supports in the flying-boat bottom. 
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The rigidity of the supports will lie somewhere "between the 
unattainable extremes of zero rigidity and infinite rigidity . The 
extreme of infinite rigidity normal to the plane of the sheet is 
one that may "be approximated in actual designs. It can he shown that 
the stress distribution in such a fixed-edge plate will, in most cases, 
be less favorable than the stress distribution in the elastic -edge 
plate. The strength of plates obtained from the theory will therefore 
be on the safe side if applied in flying-boat design. Reference might 
be made in this connection to a paper by Mesnager (reference 25) , in 
which it is shown that a rectangular plate with elastic edges of certain 
flexibility will be less highly stressed than a clamped-edge plate. This 
difference in stress may also be clearly seen by comparing the extreme - 
fiber-stress calculations by Levy (reference 19) and Way (reference 15) 
for simply supported plates and clamped plate 3 . 

The impact pressure on a flying-boat bottom in actual cases, 
however, is not even approximately uniform over a portion of the sheet 
covering several rectangular fields. Usually one rectangular panel of 
the bottom plating would resist a higher impact pressure than the 
surrounding panels, and the sheet is supported on beams of torsional 
stiffness insufficient to develop large moments along the edges. The 
high be ndin g stresses at the edges characteristic of rigidly clamped 
plates would then be absent. In order to approximate this condition, 
the plate may be assumed to be simply supported so that it is free to 
rotate about the supports. At the same time the riveted Joints prevent 
it from moving in the plane of the plate along and perpendicular to the 
supports. According to the 3 ame considerations a3 in the case of rigidly 
c lamp ed edg 03 , the result would be on the safe side. This case has 
never before been discussed and the study of such a problem seems to 
be of importance . 

For the combined bending and buckling problems 'the same consider- 
ations will hold. It is evident, however, that as soon as the side 
thrust is applied, there are displacements perpendicular to the 
supported edge 3 in the plane of the plate. Gall (reference 2 6 ) has 
found that a 3 tiffener attached to a flat sheet carrying a compressive 
load contributed approximately the same elastic support to the sheet 
as was required to give a simply supported edge (see also reference 20, 
p. 327). In combined bending and compression problems, therefore, it 
seems also important to study the ideal simply supported plates. The 
analytical expressions for these boundary conditions are formulated in 
the following discussion. 


Simply Supported Edge 

If the edge y = 0 of the plate is simply supported, the deflection w 
along this edge must be zero. At the same time this edge can rotate freely 
with respect to the x-axis j that is, there is no bending moment KL. 
along this edge. In this case, the analytical formulation of the physical 
boundary conditions is 
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Similarly, if the edge 
boundary conditions are 


< v V>o - 0 

W 2 ^X =0 


> 


J 


x - 0 of the plate is simply supported, the 


(6) 


(w) x=0 = 0 


/dfw 5 2 v \ 

Va X 2 + \2j x=0 


= 0 


Since w = 0 along y = 0, dw/dx and d^v/bx 2 must be zero 
also. The boundary conditions can therefore be written as 


(w) 


y = 0 


■ 0 




m 


= 0 


y=o 


~\ 




(7) 


Similarly, on the edge x = 0, 


( w ) x =0 = 0 



= 0 


If the plate has ideal simply supported edges, it must be free to 
move along the supported edges in the plane of the plate) that is, the 
shearing stress along the edges in the plane of the* plate is zero! 
Analytically, 



= 0 



= 0 
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or 



= 0 



= 0 


One more boundary condition is required to solve the plate problems 
uniquely, and this may “be obtained by specifying either the normal 
stresses or the displacements along the edges. 

For a plate having zero-edge compression, the normal stresses along 
the edges are zero. That is. 


or 



The strain in the median plane is 


€ 


t 

X 


du 1 / djA? 

dx + 2w/ 


(8) 



dv + 1 { cHA^ 

dy 2 'vdy/ 
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Therefore 

dtt _ c I _ 1 
dx x 2 \cbcy 

$y _ e i _ 1 ( d^Y" 
dy y 2 \byj 


and the displacement of the edgea in the x-direction is 


u = 


f 

e • -1 

J y=Constant 

x 2 Vdx/ 


dx 


while the displacement of the edgea in the y-direction is 

e » 

x=Constant L ^ 2 Vdy> 


v - 


dy 


The addition of side thrust may he expressed In terms of the 
change in displacement of the edges. 

If e x ' and €y' are expressed in terms of the stress function F, 


r 

J y=Cons 


y=Constant 


d^F 

1 dw 

fi 

. dy 2 

d 

i 

roll 



dx 


v = 


I d F &-F 1 




x=Con3tant 


ax 2 V 2 Vft 



<*y 


( 9 ) 


Clamped or Built-In Edge 


If an edge of a plate Is clamped, the deflection along this edge 
is zero, and the plane tangent to the deflected middle surface along this 
edge coincides with the Initial position of the middle plane of the plate. 



MCA TN No. 1425 


15 


If the x-eucis coincides with the clamped edge, the boundary 
conditions are 


( w ) y -0 " 0 



If the y-axia coincides with the clamped edge, the boundary 
conditions are 


( 10 ) 



If the edge is clamped rigidly against any displacement along its 
support, the strain in the median fibers must be zero along that edge. 
The boundary conditions are 





C‘X-0 ■ 

(& 

- ^ 


dyV x -o 

f dfg 


vdy 2 



( 11 ) 


The one additional condition required is again furnished by 
specifying the displacements along the edges as in equation ( 9 ) • 
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Riveted Panel with Normal Pressure Greater than That of 

Surrounding Panels 

The "boundary conditions which would approximate this situation 
are, if y = 0 is one of the edges. 



The first two expressions are those of simply supported edges, 
the third one gives the condition of zero strain along the supports, 
and the last one specifies that the displacement along the edge is zero. 


REVIEW OF PREVIOUS WORK 


The large -deflection theory of flat plates is given "by A. Foppl 
(reference 1) , and the difficulty of solving the nonlinear equations 
has "been noted "by Theodore von Kerman (reference 2). The earliest 
attempt to deal with these differential equations was, perhaps, mada 
"by H. Hencky (references 3 and 4), who devised an approximate method 
of solution for circular and square plates when the deflection is very 
large, the "bending stiffness "being then negligible . Following the 
same procedure, Kaiser (reference 5) solved the case of a simply supported 
plate with zero edge compression under lateral loading. His theoretical 
result checked closely with his experimental data. 

In the case of circular plates with large deflections, because of 
the radial symmetry, the two fundamental partial differential equations 
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which contain the linear hiharmonic differential operator and quadratic 
terms in the second derivatives can he reduced to a pair of ordinary 
nonlinear differential equations, each of the second order. For both 
the tending and the buckling problems, exact solutions are available 
(references 8 to 12) . The bending problem has been solved approximately 
by Nadai (reference 6) and Timoshenko (reference 7) and exactly by 
Way (reference 8) when the plate is under lateral pressure and edge 
monent . Wav gave a power-series solution for a rather large range of 
applied load. The buckling problem has been solved by Federhofer 
(reference 9) and Friedrichs and Stoker (references 10 to 12) . Federhofer 
cave the solution for both simply supported and clamped edges which 
yields accurate results up to values of N of about 1.25, where N is 
the ratio of the pressure applied at the edge to the lowest critic 
or Euler's pressure at which the buckling Just begins. Friedrichs 
and Stoker gave a complete solution for the simply supported circular 
plate for N up to infinity. To cover this range, they employ three 
nBthode. Each of the three methods is suitable for a particular range 
of values of N: namely, the perturbation method for low N, the 

power-series method for intennediate N, and the asymptotic solution 
for N approaching infinity. There is no solution, however, for 
the case of circular plates under combined lateral pressure and edge 
thrust. 

The exact solution for a thin, infinitely long, rectangular 
Btrip with clamped or simply supported edges was obtained by Boobnoff 
and Timoshenko (references 13 and 27), and the other cases were discussed 
bv Prescott (reference 14), Way (reference 15), Green and Southwell 
(reference 1 6 ), Levy (references 17 and 19), and Levy and Greenman 
(reference 18 ) . 

Prescott gives an approximate solution for the simply supported 
plate with no edge displacement; however, Prescott's approximation 
is rather rough. Way presented a better approximate solution for the 
clamped plates by using the Bitz energy method. Kaiser (reference 5) 
transform d the differential equations into finite -difference equations 
and solved them by the trial -and -error method. Green and Southwell 
extended the finite -difference study into finer divisions and solved 
the difference equations by means of the relaxation method. 

Levy (reference 19) gives a general solution for simply supported 
plates, and numerical solutions are given for square and rectangular 
plates with a width-span ratio of 3 to 1 under some combined lateral 
and side loading conditions. Levy and Greenman (references 17 and 10J 
extended this solution for simply supported edges to clamped edges . 

Their conditions are, however, limited to the case in which the edge 
supports are assumed to clamp the plate rigidly against rotations and 
displacements normal to the edge but to allow displacements parallel 
to the edge. They presented a numerical solution for square and 
rectangular plates with a width-span ratio of 3 to 1 under lateral 
pressure . 
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In summary the problem of rectangular plates with large deflections 
has been solved by three methods: namely, the energy method, the finite- 

-equations method, and the Fourier series method. These methods 
are briefly outlined In the following paragraphs . 


Energy Method 


The method of attack U3ed by Way (reference 15) is the RItz 
energy method. Expressions are assumed for the three displacements 
in the form of algebraic polynomials satisfying the boundary conditions; 
then, by means of minimizing the energy with respect to the coefficients, 
a system of simultaneous equations is obtained, the solution of which 
gives these coefficients. 


The energy expression for plates with large deflection is 

f(^) : 


V = 


V - qv + 6 £^2 + u x w x 2 + v y 2 + 

+ i (w x 2 + w y 2) + (y y + T y~2~ + V^r) 

+ 77 (u / + 2u y v x + v x 2 + 2u y w x w y + 2v x w x Wy)J j,dx dy (13) 


where u and v are the nondimens ional horizontal displacements and w 

pa 1 * 

is the nondime ns ional vertical displacement, q = , und the subscripts 

l6Dh 

indicate partial differentiation. In order that u, v, and w 
satisfy the boundary conditions for clamped edges. Way assumes (fig. 1) : 


u - (l - x2)(p2 _ y2)x(b 00 + b o2 y 2 + k 20 x 2 + b 22 x 2 y 2 ) 


v = (1 - x 2 )(P 2 


y 2 )y(c 00 + c 02 y 2 + c 2qX 2 + c 22 x 2 y 2) L 


(14) 


w = (1 - x) 2 (p2 . y 2)2( aoo + ^^2 + ^^2) J 

where P = u, v, w are positive In the positive directions of x 
y, z, respectively ; and ajj, bjj, Cjj are numerical constants to be 

determined later. For convenience, 1 i 3 taken to be the same as the 
power of x, and J that of y . 
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When V is minimi zed with respect to the coefficients a ij, 
t> and c.y, 11 simultaneous equations corresponding to the 
11 constants, are obtained as follows: 


o; i?~ ‘ 0 > s ?-- 0 


da. 


00 


02 


da 


20 


(15) 


dV 


dh 


00 


- 0 ’ ■ = 0 ; - ^ 2 - = 0 ; - 

' *02 *20 *22 


0 


( 16 ) 


r^L = 0; r^ 2 " = 0; r^ 2 - = 0; g^ 2 - = 0 

Sc 00 *02 S °20 *22 


(17) 


c *s 


These equations are not linear in the constants. The first three 
equations (equation ( 15 )) will contain terms of the third degree 
in the a's. Equations (l6) and (17) are linear in the b s and 
and quadratic in the a's. Way solved equations (lo) and (-*-7) for s 
and c's, respectively, in towns of a's and then substituted these 
expressions in equation (15) • There then are left three equations 
of third degree involving the a's alone. These were solved by Way 
*by successive approximations * 

Way gives the numerical solutions for cases for which fi - 1, 

1.5 and 2 for n = 0.3 up to q = 210. Since he assumed the dis- 
placements to be polynomials in x and y of finite number of terms, 
his solutions are essentially approximate. By comparing with Boobnof 
exact solution for the infinite plate, Way estimated that the error of 
his solution for 3=2 is about 10 percent on the conservative side. 
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Finite -Difference Methods of Solution 
Kaiser writes the nondimensional Von Karman equations as follows: 


V 2 S = d_ 2 w d 2 w _ . g 

\dx byj dx 2 8^2 


v 2 f = s 


: d 2 F 
cbc by dx by 


d 2 w 


Sfw c) 2 f d 2 w 

^ Sf 5 ' 8? ^ = P G 


12(1 - u 2 ) 


V 2 w = M 


vSi = 


P - P r 


"\ 




(18) 




and then transforms these five equations into finite -difference equations. 
His procedure is to assume w's at all the points and then to solve 
for S s, F's, M's, and w's. If the calculated w's do not check 
with the assumed ones, he assumes a new set of w's and repeats the 
process. The work which this involves is very tedious. In fact, as 
will P® pointed out later , when the usual method of successive 
approximations is used, the process is actually divergent. Kaiser 
the simply supported square plate with zero edge conroression 

a 4 

under a uniform lateral pressure of = H8.72. His numerical solutio: 

Eh 4 

checked with his experimental results with good accuracy. 

Southwell and Green solved four example s of the problem by means 
of a technique based on the relaxation method. The fundamental 
requirements for use of the relaxation technique are a simple finite - 
difference pattern of the variables and a simple expression of the 
boundary conditions. In using this, Southwell and Green expressed the 
( liPP eren t'ial equations in terms of the displacements u, v and w 
which then gave simple boundary conditions. Instead of using exact ’ 
relaxation patterns , they worked with the patterns which are given 
by the linear terms of the differential equations and made corrections 
from time to time, the nonlinear terras being combined with the 
"residue." 
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v\r = - + 3- 

D h 2 


du dv 1 /dv\ 2 1 / dw\ 2 t)^w dv c3u 

dx + ^dy + 2 \dx) + 2^ (jiy/ _ dy 2 [By + ^dx 

, l/^c>w\^ . 1 /dw^p . /-, \ d'w Bdu . dv . chr dw>! 

2 w + 2^s;J + (1 ■ ^ srs*^? + £ + s ^Jr 
|- (% + ^ + i-^u + ^ t- +f^1 

dx \dx dyy 1 + n 2 Sx |_v3xy \Stys 


|/|* + Ir 

dy\dx 5yy 1 + n 2 3y 


B *d) 


+ i-^JU |»yAr = 0 

1 + |i ox 

+ 1 -~- V ^V^w = 0 

1 + ja oy 


>(19) 




It is readily seen that, in order to obtain a simple expression 
for the "boundary conditions, not only is the number of the partial 
differential equations increased from two to three, but also the form 
of the terms involved becomes more complicated and the number of terms 
is Increased. This technique proves very laborious in practice. 

Equation (19) , expressing conditions of equilibrium, could have 
been derived by minimizing the total potential energy V, which is 
given by the expression 


L 2 V 

^ D = 1 + 12 + I 3 


( 20 ) 


where 




dx dy 


J 2 




+ e 


yy 


2 ^°XX e yy 



dx dy 


and 


I? = 



d-y 


where a is the lateral loading 
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The relaxation technique consists first in assuming a set of 
answers and then changing them according to the relaxation pattern 
will boundary conditions. To obtain a more rapid convergence, Southwell 
and Green multiplied the given values of w by k and substituted 
them into the energy expression to obtain 

~ | + ahl 3 (21) 


which was then minimi zed with respect to kj that is, by setting 

^ = 0 to give 
ok 


2kl 1 + kk 3 ^ - al 3 = 0 (22) 

From the third -order equation (equation (22)), k can be obtained and 
a set of values for w which are closer to the true values can be 
derived from values of k. 


Fourier Series Methods of Solution 

Levy and Greenman obtained general solutions of the rectangular 
plates (fig. 2) under combined bending and side thrust with large 
deflections by means of Fourier series. Their approach to these 
problems is given in the following discussion. 

Simply supported rectangular plates.- In order to satisfy the boundary 
conditions, w is assumed to be given by the Fourier series 


w = 


m?=l,2,3 n=l ,2,3 


. joc . Jty 
_ sin m sin n — 
m,n a d 


(23) 


The normal pressure may be expressed as a Fourier series 




r=l ,2,3 s=l,2, 3 


. joc . fly 
p sin r — ■ • sin s r~ 
^r , s a b 


(24) 
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For the compatibility equation to be satisfied 


- 2 - 2 
F = _ V *L . v iL + 

2 P = 0 , 1.2 < 1 = 0 , 1,2 


OO CO 

Y Y v 


, F must be given by 


JtX 'v / oc\ 

cos p g[- cos q 


where p* and py are constants equal to the average membrane pressure 
in the x- and y-directions, respectively, and where 


b 


P,1 



E 



(Bq + B2 + B3 + B4 + B5 + B£ + B7 + Bg + B 9 ) (26) 


and 

= Y- [ kt(p ■ *><’ ‘ ‘ ' t) 2 ]''k,t"p-k,q-t 

k=l t=l 


if q ^ 0 and p ^ 0. 


B-j^ = 0 if q = 0 or p = 0. 

B 2 = ]f ^[kt(k + p)(q - t) + k2(q - t )5 W k,t W k+p,q-t 
k=l t=l 

if q ^ 0* 


B 2 = 0 if q = 0. 


00 ol 

B 3= ^ ^>~[(k + - 4 ) + (k + p) 2 (q “ t) J w k+p,t w k,q-t 

k=l t=l 

if q ^ 0 and p ^ 0. 

B 3 =0 if q = 0 or p=0. 
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P-1 


Bk 3 zl 2L [ tt(p ' k)(t + 1) + k2 (* + <D^]w w p _^ t+q 

k=l t=l 


if p^O. 


®4 = 0 if p = 0. 


p-1 OO 

b 5 =51^1 + ^P ‘ k > + 

k=l t=l * ' 


if p ^ 0 and q ^ 0. 


B 5 = 0 if p = 0 or q = 0 . 


Be 



k=l t=l 


[kt(k + p)(t + q) - k 2 (t + <l)^v k>t v k+p 


t+q 


if q 4 0. 


= 0 if q = 0. 


00 00 

B 7 * 51 51 [**<* + ’><* + P) - ^Ck.tV'lwiv 
lc—1 t~l 


if q ^ 0 and p ^ 0. 


Bj = 0 if p = 0 or q = 0. 


OO oo 


b 8 - + r)<* + «> - < k + .»> a < t + i>dvp,t w *, 


t+q 
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If 1^0 and p / 0. 
Bq = 0 If p=0 or q = 0. 

00 00 —i 

B 9 = 2I2l[i k + p)(t + q)kt " (k + p)2t J w k+p, 

k=l t=l 


if p / 0. 


B 9 = 0 if p =» 0. 


The equilibrium equation is satisfied if 
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P r,s * (f % + » 2 s) ' Px^r.s^ ^ {1 


r b 


2^ (s - t)k - (r - *H 2 Vk,e-t w k,t 


k=l t=l 


00 00 


■ X X X* + r) • t(t + “>] 

/V X A * ' 


k=0 t=0 


00 00 

XXX + r)(t + o) • 


k=0 t=l 


- (k + r) (t + s) 


k=l t=0 


*] \+r,t w k,i 


k=l t=l 


r <x> 


(t + s)k - (k + r)tj \ +r , t+B v k>t 


tk + (r * k)(t + s) 


k=l t=0 


>J Vk.t'V 


r oo •— — 12 

X 2] L (t + a)k + {r • k)t l Vk,wk,t 


k=l t=l 


oo 8 


(a - t)(k + r) + tk 


k=0 t=l 


*] ^k,e-t w k+r,- 


oo g 

+ X XX • t)k + *<* + r )]v,»*t 

k=l t=l ' ' 
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When the lateral pressure is given, Pp 6 can he determined. 

Equation (27) represents a doubly infinite family of equations. In 
each of the equations of the family, the coefficients b p ^ may he 
replaced hy their values as given by equation (2o) . The resulting 
equations -will involve the known normal pressure coefficients p^g, 
the cubes of the deflection coefficients n , and_the average 
membrane pressures in the x- and the y — directions, P x and Py > 
respectively. Values of p x and p y can be determined from the 
conditions that the plates are either subjected to known edge compressions 
or known edge displacements. The number of these equations is equal 
to the number of unknown deflection coefficients w rrl ^ n . 

The procedure now is, with the known values of p r ^ a , P° assume w p p 
and to solve the other coefficients by successive approximation. However, 
the work involved is tremendous, and it is very easy to make mistakes. 

As illustrated by Levy in a relatively simple case of a square plate, 
if six deflection coefficients are used, then each equation contains 
60 third -order terms. And for each given applied normal pressure 
these six 60-term, third-order equations must bo solved by successive 
approximations . 

Clamped rectangular plates.- Levy and Greonman solved the case of the 
clamped rectangular plate by assuming that the edges are clamped rigidly 
against rotations and displacements normal to the edges but are permitted 
to move freely parallel to the edges . 

The required edge moments m^ and m^ are replaced by an 
auxiliary pressure distribution Pa(x,y) near the edges of the plate. 

The auxiliary pressure can bo expressed as a Fourier series as follows. 


Pa( x >y) 


I" 

r =1,3,5 a ' 


TTOC 

r sin — - + 


z 

8^1,3, 5 


4jtrn x 

b2 


s sin 


532 

b 


(28) 


By writing and nu as Fourier series, where 

coefficients to be determined, 


“x 



OO 


r =l,3,5 


kr 


sin 


ntx 


^s 


and kp 


are 


"V 




s-l, 3 ) 5 



> 




(29) 



28 


MCA TN No. 1425 


Inserting equation (29) into equation ( 28 ) gives 

P 00 00 

Pa(*,y) -Q) P X ( rk s + a**-) sin sin £2 ( 30 ) 

r-1,3,5 s=l, 3, 5 


On combining the auxiliary pressure p a (x,y) with the normal 
pressure p z , equation (24), the following equation is obtained: 

OO 00 

Pc( x ,y) = 23 23 sin r ~ sin s ^ ( 31 ) 

*=1,2,3 s=l,2,3 


where 



(rpk 8 + spkp) + 


(32) 


Since the edge moments m^ and n^- have been replaced by the 
auxiliary pressure distribution p a (x,y), the general solution for the 
simply supported rectangular plate (equations (23) to ( 27 ) ) can be 
applied to clanqjed plates, and the remaining problem is to determine 
the values of k g and kp. These values are obtained by use of the 

boundary condition that the slope at the edges of the plate is zero. 
Equating to zero the normal slopes along the edges gives 



00 



nt=l,3,5 



*=1,3,5 


f cHm 

vty/y=0,y=b 


» 0 


CO 


r. 

m=l,3,5 


00 


X 

n= l>3,5 


22 

b 


m,n 


, nuac 
sin 

a 


> 


( 33 ) 


y 
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Equation (33) is equivalent to the set of equations 


A 


0 = w 1;1 + 3 w 1<3 + 5w 1>5 + 


° * "3,1 + 3w 3,3 + ^3,5 + 


> 


0 = w_ , + 3w s o + 5W- c + • • • 

0 y*- y y o y ) y 


j 


(34) 


The deflection coefficients w m n must now "be solved from the 
family of equations (equation (27)) £or the linear term in terms of the 
cubic terms and the pressure coefficients ? r s . The express long for v m ^ n 
thus obtained are now substituted into equation (34), and the expressions 
for pressure coefficients p r s are obtained from equation (32). The 
resulting family of equations’ contains linear terms of pkp and pk g 

and. the cubes of the deflection functions w m ^ n . 


The method of obtaining the required values of the deflection 
coefficients w m n and the ed©3 -moment coefficients pk^. and p^ g 
* Wf -j_ pa 4 

consists in assuming values for - 7* J and tnen .solving for , , 

h Eh 4 


pk g> pkp, ... by successive approximations from the 


v l,3 
h 7 * 

simultaneous equations. The procedure is even more laborious than that 
for 3 imply supported plate 3 . Two numerical solutions are given, namely 
solutions of the bending problem for a square plate and for a rectangular 
plate with length-width ratio of 1.5. 


FINITE -DIFFKRSNCdl -SQUATIOMS OF BOUNDARY -VALU3 PROBLEMS 


Some fundamental concepts about the f toite “difference approximation 
may be worthy of mention before the partial differential equations are 
converted into finite -difference expressions. 

It i 3 assumed that a function f(x) of the variable x is defined 
for equidistant values of x. If x is one of the values for which f(x) 
is defined, f (x) is also defined for the values of x + k Ax, where Ax 
is the interval between two successive values of x and k is an integer. 
For the sake of simplicity, the value of the function y =f(x) for 
x + k Ax my be written as: 


f(x + k Ax) - ^x+kAx 
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The first difference or the difference of the first order Ay x 
of y at the point x is now defined as the increment of the value 
of y obtained in going from x to x + Ax: 


^x ^x+Ax ” ^x 


It is seen that the increment in the direction of increasing x has 
been arbitrarily chosen; Zty x could also be defined by the difference 
y x - y^-Ax • This process is continued and the increment of the first 
difference obtained in going from x to x + Ax is called the 
difference of second order of y at x; that is, 

A ^x = A ^X4 Ax “ ^x 


(V x+2Ax “ ^x+Ax^) 



. - yx+2Ax “ ^x+Ax + y x 


In general, the difference of order n is defined by 

A n n-1 A n-1 

A = A yx-tAx - A y x 

If Ax is chosen equal to unity, 

y x+nAx " y x+n 


By the use of this notation, the sequence of differences becomes 
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- y x+ l ’ y 3 


A2y x y x+2 ' 2y x+l + y x 


A y x y x+3 ~ 3y x + 2 + 3y x+ l y x 




r=0 


n: 


r!(n - r) I 


'x+n-r 


(35) 


In many physical problems only differences of even order occur. 

In such cases It is more convenient to define the differences 
2m 

A y in the following way: 

A y x ^x-l ~ ^x + y x+l 

p 

That is, A y x is the increment of the first difference taken on the 
right- and left-hand sides of the point x. In general, 

A 2m y x = (36) 

In this case a difference of order 2m represents a linear expression 
in yx-nu yx-nvfl> y x * y x+m-l> y x-nn* 

In replacing partial derivatives by the finite -difference 
expressions, the differences corresponding to the changes of both the 
coordinates x and y are considered. With the notations as shown 
in figure 3, the first differences at a point Pl^ n in the x- and 
the y-directions are, respectively: 


A r w m,n w m+l,n ~ W m,n 
‘ ! V w m,n = w m,n+l " w m,n 
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The three kinds of second differences are as follows: 


A_.w = A v 
m.n jc m,n 


w m + l,n - *k,n + w m-l,n 


/ ^yy w m,n ^y W m,n 


v - 2w_ _ + v _ . 

m,n+l m^n-l 


W = Vm+l,n - Vm,n 


(^m+l,n+l w m+l,n) (ym,n+l w m,i^) 


= V “W - v +v 

m+l,n+l m+l,n m,n+l m,n 


The three kinds of fourth differences, which will he used later, are : 


\m w m,n \ W m,n 


w „ - 4w + 6w 

m+2,n m+l,n m,n 


^ v m-l,n + w m-2,n 


&yyyy v m,n “ ^y w m,n 


= w ™ - + 6w_ _ - 4w_ _ i + w„ „ 0 

m,n+2 m,n+l m,n ni^n-2 


A w = A _ 

xxyy m,n xy m,n 


= v„,i - 2w . + v , , - 2w , + 4v 

m+ l> n+ l m+l,n m+l,n-l m,n+l m,n 


- w _ , + w - 2 w„ i + w -i _ -i 

m,n-l m-l,n+l m-l,n m-l,n-l 
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Partial derivatives may "be approximated "by finite differences as 
follows : 

dw ^ ^ dw _ ^y v 

c*x Ax * Sy 

2tx 2 Ax 2 * C»y 2 Ay 2 

_£sl a ^sHL > (39) 

5y Ax Ay 

d\, A . V* 

dx 1 * oy^ A? 

jv _ V 2 * 1 

3x 2 dy 2 Ax 2 Ay 2 
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If Ar = £y = Al, and. the relations (37) and. ( 38 ) are used, 
equation (40) may be written as 


F m+2,n ®m+l,n + 2 ® F m,n ” ®m-l,n + F m-2,n + F m,n+2 “ ^n^n+l 


- ■, + F„ 


,n-l + m,n-2 + ^uH-l^n+l + ^m+l^n-l + ^m-ljn+l + ^m-l^n-l 

(^m+l,n+l " w m+l,n " w m,n+l + v m,n) ~ (*m+l,n " ^m^n + w m-l,n) 
x (^m,n+l " ^m^n + w m,n-l) 


and 


w m+2,n ” ^m+l^n + 20w m,n " ^m-l,n + w m-2,n + w m,n+2 “ ^n^n+l 

^m^-l + w m,n-2 + ^m+l^n+l + ^m+l^n-l + ^m-l^n+l + ^m-l^n-l 


= 10.8(AI) p + 10.8 


(*m,n+l " ^m^n + F m,n-l) 


x (*m+l,n " ^m^n + w m-l,n) + (* 111 + 1 , n ” + F m-l,n) 

X (^m,n+l ” ^m^n + w m,n-l) " 2 (*L+ 1 , 


n +1 ~ F m+l,n " F m,n +1 + F m 




x (>iy 


n +1 w m+l,n w m,n+l + 


- Yl 

m,nj 


(^ 2 ) 


In actually writing these equations for each net point, it is more 
convenient to employ the finite -difference pattern or so-called 
relaxation pattern as shown in figure 4 rather than to substitute 
directly into equations (4l) and (42) . 

In terms of finite differences, the boundary conditions can be 
formulated in the manner discussed in the following paragraphs. 
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Simply Supported Edge 

The 'bo undar y conditions for th© simply supported edge y = 0 are : 



and, for plate b with, zero edge compression: 



or, for plates with zero or known edge displacements: 


iL 




to? ^ 





dT = v 


Let n = 0 denote the edge points along y = 0. The finite - 
difference expressions for the "boundary conditions are : 





J 


36 


MCA TN No. 1425 


for plates with zero edge compression and 


Z [v* - ^ - 1 (v)‘ 


= v. 


li.n 


where n = 0 and n = k denote points along the two edges y = 0 
y = "b, respectively, and i denotes any point along the line z 
in the plate . 


Clamped Edge 

The "boundary conditions for the clamped edge y = 0 oj^e : 


(v) 




ar 0 


C 6 2 f \ 
V* 2 ' 


= 0 


n 



u 



l/'dwN 2 

V ‘ 5 W 


dy = v 


With the same notations as were used for the simply supported 
the finite -difference expressions are: 


(45) 


and 

= Constant 


edges. 
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Riveted Panel with Normal Pressure Greater than That of 

Surrounding Panels* 

The "boundary conditions which approximate this case are • 


(v) 


j=0 


0 



if y = 0 is one of the edges. 

Expressed in terms of finite differences f these conditions 
"become : 
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v m,0 * 0 | 

($%,o ■ 0 

- 0 

k-1 r 

51 " |(V> 2 “ 0 

n=0 * i i,n ^ 

t 

The boundary -value problem which approximates the riveted sheet- 
stringer panel subjected to uniform normal pressure higher than that 
of the surrounding panels may be formulated in terms of finite 
differences. 

In order to start with a simpler case, the square flat plate will 
be discussed, since, on account of symmetry, only one -eighth of the 
plate need be studied. 

The finite -difference approximation of any differential equation 
requires that every point in the domain to which the equation applies 
must satisfy the Initial differential equation. If the points to be 
taken are infinite In number, the solution of the difference equations 
is the exact solution of the corresponding differential equations . But 
the points to be taken are finite In number, the solution will be 
approximate, and the degree of approximation will Increase as the 
number of points taken is reduced. 

Since the diagonals of a square plate are axes of symmetry. If 
the boundary conditions along the four sides are the same, k = v k ^ 

and € i] c = e ] £ i» The conditions for zero edge displacements may 
be put Into different forms. Since 





In terms of finite differences 
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€ x + f y 


2 2 F _ d% + 0% _ 0% 

dy 2 cbc 2 ox 2 dy 2 


■(9 * 3> • ■’ 


= [^ 2 F) (1 - n) 


Note that 


v 2 = + 


2k 2 §y 2 


Equation (48) then "becomes 


i,0 + 2 ( v2p ’)i,l + **• + 2 ^ 2 *)i,n + •“ + 2 ^ 2 ^i, k -l 


iA (1 - h) 


v i,n+l " w i,n 


This simplification is not necessary, "but it is useful in applying the 
relaxation method. 

n__=_l.- On referring to figure 5, it ia seer, that points 1' and 2* 

are fictitious points placed outside the plate in order to give a 
"better approximation to the "boundary conditions. 

By the use of |i 2 = 0.1 or n = O .316228 for aluminum alloy, 
the compatibility equation "becomes 


r 0 - 32F-L + 8 f 2 + 4f 1 = Kq 


where 


• P 2 

K 0 = ( w 2 " 2*1 + v 0 ) - ( 2 w x - 2 w 0 ) 
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Then the equilibrium equation is 


20w 0 - 32W! + 0w 2 + 4w 1 , = p' + 10.8^2(2F 1 - 2F Q )(2w 1 - 2w Q ) 

- 2(w 2 - 2w ± + w q )(F 2 - 2F X + F 0 ) 

where p' = 12(l - n 2 )(Al)S = 0.6?5p, since AZ = 

The boundary conditions are : 


(51) 


(a) = 0, w 2 = 0 

(h) w lt - 2w x + w 0 = 0 

(c) F q - 2F X + F-j_i - h(2F 2 - 2F X ) = 0 


(d) (4F X - 4F q ) + (F 0 + 2F 2 + F 1 , - 4F X ) = S x 


t 

where = (w-^ - Wq)‘ 


(vi - w 0 ) 2 


• The "boundary -value problem 


(1 - n) 0.341886 

now determines the values of w uniquely and the values of F to within 
an unknown constant. Since the actual value of the constant is irrelevant, 
it may he defined by letting F 2 = 0. 

Solving w., , , w 2 . and F,i from the boundary conditions gives 

the following result: 


w l' = -^o 

w 2 , = -w-l = 0 

F l* = ~F 0 + 2(1 - \x)F ± 
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When these values are substituted into equations ( 50 ) , (51) > and (d), 
the resulting equations are 


16F q - 26.529824F 1 = -3 w q 2 


16 w q = p' + 43 • 2 w q F 0 


-4f q + 1.3675^^ = 




0.341886 


(52) 


The eight or nine significant figures in these equations are due to 
computations made with a computing machine having 10 columns. In order 
to get satisfactory results in subsequent computations it is convenient 
to retain a number of figures beyond those normally considered justi- 
fiable because of the precision of the basic data. 

n = 2 . - With reference to figure 6, points 3' > V , and 5’ are again 
fictitious points. The compatibility equations are: 


2 OF - 32F + 8 F + kF = K 

0 1 2 3 0 


- 8 f 0 + 25 F l - i 6 f 2 - 8 f 3 + 6 f u + F , = K x 


2F - 16F + 22F + 4F - 16F, + 2F + 2F. . = K 

0 1 23 45 4’ 2 


> (53) 


(\y v ) 2 - A x 2w *y 2v 


where Kq, and are equal to 

1 , and 2 , respectively. 


at points 
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The equilibrium conditions are: 


20vq - 32w^ + 8w 2 + 4w^ 

= p' + 21.6 JJoq' + Pq')^ - V 0 ) - r'(w 0 - 2w x + W 2 )J 
-8 wq + 25w^ - lowg - 8W3 + 6 wi 4 . + w^i 

= p’ + 10.8 jo.i'(2w 2 - 2w^ + Pi'( w o " 2Wi + v ^) 

" 2V K " w 3 - v 2 + w x )J 
2 wq - l6vj_ + 22 w 2 + 4w^ - l6v^ + 2w^ + 2wj + , 

= p ' + 10 . 8 (a^ • + pjj) (w^ - 2 w 2 + w^ - 27 2 ' (wtj - 2w^ + Wg^ 


where a', 3', y 1 are Z^F, A 2 F, A^^F at the respective points 

indicated by the subscripts. 

The conditions for zero edge displacements are: 


-2Fq - 3Fi + 4F 2 - 2 F 3 + 2Fjj. + F 1 = S-l 


F 0 “ 5F 2 + 2F^ + F^ + F^i = S 2 


where 


1 0.341886 _f Wl w °^ + ^ -1 

3 2 * 0-341886 [ (W 2 ' W l )£ + K ' W£)2 ] 
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boundary conditions are: 



(a) 

w 3 

II 

O 

£ 

II 

o 

N* 

V 5 a 

0 

0») 

V 

- 2w 3 + v-]_ = 

0 



w 4 , 

- 2w 4 + w 2 = 

0 



w 5‘ 

- 2vr 5 + W 4 = 

0 


(c) 

F i 

- 2F 3 + F , - 

d(2F 4 

- 2F 3 ) = 0 


F 2 - 2F 4 + F v - h(F 5 - F 4 + F 3 ) = 0 
F 4 - 2F^ + F,_, =0 

For the same reason as explained in the case of n = 1, let F^ = 0. 
Solution of the houndary -conditions equations gives 

(d) w^» = -w^ 

V = "*2 

w 5 , = 0 

(e) F 5 - = -F 4 

F 4 , = 2F 4 + h(F 3 - 2F 4 ) - F 2 
Fy = 2F 3 + h(2F 4 - 2F 3 ) - F 1 
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The combination of the foregoing equations (d) and (e) with equations (53), 
(54), and (55) gives 


20Fq ~ 32F^ + 8F 2 + 4F ^ 13 Kq 

-8f q + 24f x - i6f 2 - 6.632456f 3 + 6.632456F4 - % 

2F 0 - 16f x + 20F 2 + 4.632456f 3 - 13«264912F^ = Kg j> 
-2Fq - 4Fjl + 4f 2 - 0.632456f 3 + 2.632456f^ - S x 
F 0 - 6 f 2 + 2 . 316228 f 3 + 1.367544F^ « Sg 

J 


(56) 


[20 + 21.6(0^ + P 0 * + 7 q ' )j wq 

- [32 + 2l.6(a 0 ’ + 3 0 * + 27q ' )j w x + 

“(8 + XOeSP-j^’JVQ + + 2X»6(ctr^ 1 + 3^ 1 

- [16 + 12. 6(0^ + 7 1 ')jw 2 = p' 

2w 0 - jl6 + 10.8(02' + 3 2 ^] w 1 

+ (20 + 21.6(02' + 3 2 ' + 7 2 ')jw 2 = p* 

J 

where p* = 12(l - n 2 )( Al)\ = 0.042l875p, since A1 = i. 

n = 3 * - Reference is made to figure 7 and to the fact that points 6', 

7', 8', and 9' are fictitious points for reasons explained in the case 
n a 1 ^ then the compatibility equations are as follows : 


(8 + 21.67o')w 2 = p' 
+ 7 1 , )]w 1 


> (57) 
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20F 0 - 32F l + 8 F 2 + 4F 3 = Kq 

-8f 0 + 25*1 - i6f 2 - 8f 3 + 6f 4 + Fg = K-l 

2F q - l 6 F 1 + 22 F 2 + 4F 3 - 1oF 4 + 2 F^ + 2 F ? = 

F 0 - Sf-l + 4f 2 + 20F 3 - i6f 4 + 2F 5 - 8f 6 + 4f t + f 6 , = ^ 

3F-,_ - aF 2 - QF 3 + 23F 4 - 8F 5 + 2F 6 - 8F 7 + 3Fg + F ? . = K 4 

2F 2 + 2F 3 - l6F 4 + 20F^ + 4F 7 - 16 Fs + 2F^ + 2Fg. = K^ 



where Kq, K x , Kg, K^, and are equal to 

at points 0, 1, 2, 3, 4-, and. 5, respectively. 
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The equilibrium equations are : 


20wq - 32w^ + 8w 2 + 4wg = p' + 21.6 ^(oo' + P 0 ')( w l " w 0) 

- V< w 0 ' 2*1 + w 2^] 

-8w 0 + 25^1 - l6w 2 - Qw^ + 6w 4 

= p' + 10.8 - w-j^) + P 1 '(w 0 - 2w x + Wg) 

" 27 1 , (w 1 - v 2 - w 3 + w 4 )J 
2w q - 16 W-L + 22w 2 + 4vg - 16 w 4 + 2v^ 


P ' + 10 ‘ 8 [(<*2 ' + V * < W 1 “ 2 v 2 + w 4 ) - 27 2 *( w 2 


2W], + 


= p' +10.8 a 3 , (2w 4 - 2 w 3 ) + P 3 , (v 1 - 2w 3 + wg) 

" 27 3 '( w 3 - W4 - w 6 + v ? )^ 

3»i - 8 w 2 - 8*3 + 23v 4 - 8w 5 + 2wg - 8w ? + 3 w 8 + w ? . 

= p' + 10.8 a 4 '(w 3 - 2W 4 + w 5 ) + P4'(w 2 - 2W4 + 

~ 27 4 '(w 4 - w 5 - w ? + wg)] 

2w 2 + 2w 3 - l6w 4 + 20 w 5 + 4 w - l6wg + 2w^ + 2wg. 

= p* + 10.8 (05' + p 5 .')(w 4 - 2w 5 + wg) - 27 5 '(w 5 - 2wg + V9) 
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where a', 0*, and y' are A^F, A^F, A^F at the respective points 

corresponding to the subscripts, and p' = 12(1 - n 2 )(Al)^p = 0.00833333p, 

since Al = 7. 

o 

The conditions for zero edge displacements are : 


- 2 F 0 - 2 F X + 4F 2 - 5F 3 + 4F 4 - 2F 6 + 2 F 7 + Fg. = S 2 

F 0 - 4F 2 + 3F 3 - 3F 4 + 2F 5 + Fg - 2F ? + F q + f t , = S 2 

F x + 2F 2 - 2 F 3 - 2 F^ - 5F 5 + Fg + 3F 7 + F^ + Fq. = S 3 








where 


(60) 


s * - rb 5 k+v 


The "boundary conditions are: 



(a) v 6 =0, w ? = 0, v 8 = 0, w 9 = 0 

(b) wg, - 2wg + w^ = 0 

w^, - 2w^ + = 0 

wg, - 2wg + w^ = 0 


- 2w^ + wq = 0 
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(c) F 3 - 2F 6 + F 6 , - n(2F ? - 2F 6 ) = 0 

F^ - 2F-j + F^, - u(Fg + Fg - 2F-j) =* 0 

F 5 - 2F 8 + Fg, - h(F 9 + F ? - 2F q ) = 0 

V - 2F 9 + Fg = 0 

Solutions of the boundary -conditions equations give 

( d) w 6 , = -w 3 

w 7 ' = " w 4 

v 8* = " w 5 

v = 0 

9’ 

(e) Fg, = -Fg + 1.367544 f 6 + 0.632456f ? 

F , = -F^ + 1.367544f 7 + 0.31622^ + 0.31622^3 

F gt = -F 5 + 1.36?544 f 8 + 0.3lS228F T 

V - " F 8 

where = 0 is assumed for the same reason as explained in the case 
of n = 1. 
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Combination of the foregoing equations gives: 


r ! + ^2 + ^3 


-&0 + 25Fi ' i6f 2 * ^3 + + f 6 a k 1 


2F q - 16F! + 22Fg + - l6F^ + 2F^ + 2F^ = 

F 0 - 8Fi + 4 f 2 + 1^2 " I0F4 + 2*5 - 6.6 32456F5 + 4.632456F7 = 
3®i ’ ^2 " ^3 + 22F 4 " ^5 + 2.31622^6 - 6.632456 f 7 
+ 3*31622^8 = K4 

2 F 2 + 2F 3 - loF^ + 1^ 5 + 4 . 632456 f 7 - 13 . 264912 F 0 = &- 
- 2 F 0 - 2 Fi + 4 F 2 - 6f 3 + 4 F^ - O.632456F6 + 2 . 632456 f 7 =» S 1 

F 0 - 4 f 2 + 3F3 - 4 f^ + 2F5 + I.316228F6 - 0 . 632456 f^ 

+ 1 . 31622^8 a S 2 

F i + 2F 2 - 2F3 - 2F4 - 6 f 5 + Fg + 3*3 i 6228 f 7 + 1.367544FQ = s 3 
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and. 


» 3 ~ |j t ' j + 10.8(a 2 ' + P 2 ' + 47 2 *)Jv4 


J20 + 21.6(0^' + p 0 * + 7o')jwo " j32 + 21.6(00* + Po* + 7 0 1 ) 

+ ^8 + 21 . 67 0 ^v 2 + = p' 

- (^8 + IO.SP^wq + 25 + 21.6(0^ + P^ + 7^*) v-^ 

- jlo + 21.6(0^' + 7 X * ) w 2 - ^8 + 10.8^* + 27 1 , )jw 3 

+ (6 + 21.67 1 ')w 1 ^ = p* 

2w q - [l6 + 10.8(a 2 ' + Pg’^l + [22 + 21.6(02* + p 2 * + 7 2 ')]w 2 
+ 4w 3 - |l6 
+ (2 + 21 . 67 2 ')w 5 = p’ 

v Q - (j3 + 10 . 8p 3 J^ w x + 4 w 2 + jjL9 + 21.6(03 ' + P3' + 7 3 ')Jv 3 
- jjLo + 21 . 6(03' + 73')] v 4 + 2 w 5 = p' 

3W-L - (8 + 10.8P4')w 2 - (8 + 10.8a^')w 3 

+ 22 + 21.6(04* + p 4 * + 74^4 - |8 + 10.8(04* + 27 4 ') 

2v* 2 + 2W3 - fl6 + 10.8(0^' + P5*)jw4 

+ pL8 + 21.6(0^' + P 5 * + 7 5 ')Jw 5 = p' 


V(62) 


w 5 = p 


J 
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METHOD OF SUCCESSIVE APPROXIMATIONS 
Explanation 


After the "boundary -value problems are expressed in terms of 
finite -difference equations, two sets of simultaneous equations are 
obtained. The first set consists of the compatibility equations the 
equations specifying the condition of zero edge displacements. These 
equations contain linear terms of the nondiraensional stress function F 
and the second-order terms of the nondiraensional deflection w, and 
are of the form 


c F + c F + ... + c^ F = K-. 
00 0 01 1 On n u 

C 10 F 0 + C 11 F 1 + “• + c ln F n = IC L 


and 


c Vo + C 'll F l + ••• +c, m F n = s i 


J 


where K 1 = (Axyw) 2 - (a^-w) (A y 2 w) at points 0,1, and so forth, 
corresponding to the subscripts of K; = — - — 

i - H A_ 

and- c oo> c oi> •••> c 'io> c 'n> ••• are given constants. 


( 63 ) 


The second set consists of the equilibrium equations, which contain 
the linear terms of w with coefficients involving linear terms in F 
and are of the form 


(*00 + ^OO^O + fc'ooP'o + * "oo 7 ’o) w 0 

+ (foi + *01*0 + ^'oiP'o + b Wo) w l 
+ ... + + ■b 0 n a 'o + ^'onP'o + ^Wo)^ = P* 


( 64 ) 
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where a’ = A^F , P' = Ay 2 F, 7 1 = at points 0,1, • • • 

corresponding to the subscripts of a', 3', 7', and Uqq, sqi, • • •, 
boo> b>oi, • • •, b ’qq, "b , 01 , . . ., b"QQ, b 0 i> • • ♦ are g iven 
constants . 

If a set of values of w is assumed at each of the net points 
and the values of K ± and Si are computed, equation (63) becomes 

a system. of linear simultaneous equations in F and can therefore 
be solved exactly by Crout's method for solving Bystems of linear 
simultaneous equations (reference 24). After the values of F have 
been computed from equation (63), values of a', P 1 , anj i ~f x can 

be found" without any difficulty. Then equation (64) becomes another 
system of linear simultaneous equations and may be solved exactly by 
Crout's method again. If the values of w found from equation (64) 
check with those assumed, the problem Is completely solved. 

In most cases, however, the values of w will not check with 
each other. By following the usual method of successive approximations, 
the computed w's will now replace the assumed ones and the cycle of 
computations will bo repeated. If the value of w at the end of the 
cycle still does not check with the one assumed at the beginning of the 
cycle, another cycle will be performed. In this problem, however, if 
the ordinary method were followed, the results would be found to diverge, 
oscillating to infinity. Therefore, a special procedure must be devised 
to make the process converge. 

A simple case will be examined first. In the boundary -value 
problem in which n = 1 under the normal pressure p = 100, equation (52) 
can easily be reduced, to the form 


w 0 = 


16 + 37.6903 w 0 2 


or 


w 0 3 + 0. 42450 7w q - 1.790888 = 0 


(65) 


The third-order algebraic equation can easily be solved, and the roots 
of this equation are 


w Q = 1.098254 and (-0.549127 ± 1.1528781) 
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For the physical problem, only the real root is of interest because 
the imaginary roots do not have any physical meaning. 

An attempt will now be made to solve equation (65) by the usual 
method of successive approximation*. It is assumed that 


w 0 

w 0 


1.200000, w 0 2 = i.44oooo 

■■ = 0.960516 
70.27474 

0.922591 


If it is assumed that wq 2 = 0.922591 for the second cycle and that 

the value of w^ 2 found from the second cycle is the value for the 

L p 

third cycle, and so on, the following values of Wq are found from 

various cycles: 


1.767416, 0.667554, 2.689324, and so forth. 


These values are oscillatorily divergent. A plot of these values 
against cycles shows that they oscillate about the true value 1.206l6l, 
and the true value Is approximately the mean of the values obtained 
from two consecutive cycles (fig. 8). 

If w 0 2 = Jf 1 - 1 ^ 000 + 0.922591) = 1.181296 is taken as the 

assumed value of wq 2 for the second cycle, and the mean of this 

value and the value found from the second cycle are taken as the 
assumed value for the third cycle, and so forth, the values of Wq 2 

are found from various cycles as follows: 



NA.CA TN No. 1425 


55 


Cycles 

2 


3 

4 

5 

p 

Wq assumed 

1.181296 

1.212550 

1.204658 

1.206524 

2 

Wq found 

1.243805 

1.196766 

1.208390 

1.204526 

Cycles 


6 


7 


Wq^ assumed 

1.206075 

1.206182 


w q 2 found 

1.206289 

1 .206131 



This process is convergent and w Q converges to the real root of 
equation (65) • The value of w 0 found at the end of the seventh cycle 
is 1.09824-0 and is accurate to four figures at the end of the fifth cycle, 
in which case it is found to he I.O98OIO. The results are plotted 
against cycles in figure 9* 

It is to he noted that Kq = -3w Q 2 in the case of n = 1. The 
values obtained hy the method of successive approximations would 
converge if Kq were assumed to he the mean value of two consecutive 
cycles. It is found that this convergent property is the same for n > 1. 
If the mean of K's or S's found from two consecutive cycles is 
taken, the values are convergent hut are oscillatorily divergent if the 
usual way of successive approximations is followed . 

It may he pointed out here that for the special case n = 1, if 
the mean of the values of w^ from two consecutive cycles is used, the 
values are also convergent, and if Wq 2 for the second cycle is 
assumed to he equal to the sum of 0.6 times the assumed value for the 
first cycle and 0.4 times the value found from the first cycle, and 
so on, the convergence is much more rapid (fig. 10) , hut this result 
is not true for the cases with n > 1. 


The rapidity of the convergence depends on the accuracy of the 
assuned values of K's and S's for the first trial. The deflection w 
from the linear small -deflection theory can easily he determined. When p 
is small, the values of w so determined would give a first approximation 
to the problem. It is convenient, therefore, to start the computation 
when p is smal 3 and then to consider the cases when p is large . 

Also it is advisable to begin with hut a few net points and then 
gradually to increase the number of net points. For example, consider 

case n = 1. When w 0 is found for a certain small p, a 

curve of w Q against p can he plotted because the slope of the curve 

at the origin can he determined from the small -deflection theory. For 
a larger value of p, wq can now he estimated hy extrapolation. For 
n = 2, the value of vq found for n = 1 can he used as a first trial . 
However, Wp and w-^ are still difficult to estimate. In order to 


obtain first approximations to these quantities, the ratios Wp/w ( 


and 


W 3 /W 0 


2 ' 0 

may he found from the small -deflection theory and the values 
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of w 2 and w^ computed "by multiplying these ratios "by the estimated 

value of vq. When the deflections have been assumed at every point of 
the net, the values of K and S can he computed . These are the values 
which may he used as a first trial. By successive approximations, the 
true values of the w's are then determined. The values of and 

the (v n /v Q ) 's are now plotted against p to estimate the corresponding 

values at a larger p. The values estimated hy extrapolation may he 
used as the trial values corresponding to that p. The process is 
repeated until the maximum p is reached. For n = 3, w Q from n = 2 

is used as a first trial) the remainder of the procedure is the same 
as before. 


Sample Calculations 

Finite -difference solutions of small -deflection theory .- The small- 
deflection theory of the simply supported square plate will he studied 
first. The differential equation is 

A. | (66) 


and the boundary conditions are 


w 

= 0 

along 

0 w 


along 

r\ 

= 0 

cbc 2 


0 w 

x 2 

= 0 

along 


= +* 




( 67 ) 


J 


where a is the length of the sides. 


With equations (66) and (67) written nondinensionally hy letting 
4 


W = h> 


pa 
“T> 
Sh 4 


x 1 = — . and 


= r-, where w 1 , 


and 


and y* are nondimensional deflection, pressure, and lengths, respectively, 
and with the primes dropped, the boundary- value problem is: 
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V*v 

= 12(1 - 

i-i 2 )p 

V 

= 0 

at 

X = 

"2> 

cfw 

= 0 

at 

x = 4 

<bc 2 


r\ 

= 0 

at 

y = ^ 




By retaining the notations previously used, the finite -difference 
equations for the problem are 


A^v + a^ 2 * + A^v = p’ 



where p* = 12(l - p 2 ) ( Al ) ^p . 

For n » 1 (fig. 5), the finite -difference equation, after the boundary 
conditions are employed, becomes 


l6w r 


= P 1 


theref ore , 


Wq = 0.0625P 1 
= 0.042l86p 


for p 2 = 0.1. 
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For n = 2 (fig. 6), the finite -difference equations, after the 
■boundary conditions are inserted, "become 


20 wq - 32*i + d*2 = P' 

-8wq + 2k*i - l6w 2 = p' S>- (69) 

?Wq - 16*1 + 20w 2 - p' 

.J 

When Crout's nethod is used to solve these equations the solutions 
of equation ( 69 ) are 


where 


v Q = 1.03l250p' 
w 1 = 0 . 750000p * 

w 2 = 0.546875 p' 

is takBn to he equal to 0.1. 


= 0.043506p 
= 0.03l64lp 
= 0.023071p 
For n = 0.3, 


w Q = 0.032989P 


For n = 3 (fig* 7 ), the finite -difference equations, after the 
houndary conditions are employed, "become 


20w q - 32w 1 + d* 2 + 4w^ = p* 

-8wq + 25w 1 - l6w 2 - 8w^ + 6W4 = P* 

2 *q - l6w 1 + 22 w 2 + 4w^ - l6wj^ + 2w^ = p' 
Wq - 8w 1 + 4 w 2 + 19w^ - lSwj^ + 2w^ = p* r 

3Wp - 8 w 2 - 8w^ + 22w^ - a p* 

2*2 + 2w^ - l6w^ + l8w^ = p* 

) 


( 70 ) 
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The solutions of equation (70) are: 

Vq = 5.2U6672p' 
w-l = 4.597633P ' 
w 2 = 4.031250p' 
w 3 = 2.735207P' 
vj^ = 2. 40236 7p' 
w 5 = 1. 439l64p * 

if is assumed to he 0.1. If |i 


= 0.043722p 
= 0.038314? 

= 0.033594p 
= 0.022793? 

= 0 . 020020 ? 

= 0.011993? 

is assumed to he 0.3, the answer is 


Vq = 0.044208? 

Timoshenko gives the exact value of v Q for a simply supported 
square plate (reference 27) as: 


Vq = 0.0443p 


Therefore the solution "by finite differences vith n = 3 I 0 Q rror 
hy 0.23 percent. This solution is seen to he sufficiently accurate for 
engineering purposes . The agreement of the finite -difference approxi- 
mation vith the mere exact results of Timoshenko is sufficiently close 
to encourage application of the finite -difference approximation to the 
problems vith large deflections. 

The large -deflections problem, n = 2 .- After the boundary conditions 
are inserted, the two sets of finite -difference equations are: 


20F 0 - 32Fq + 8F 2 + 4F 3 = Kq 


-8Fq + 24f x - l6F 2 - 6.632456f 3 + 6.632456f^ = 
2Fq - 16FJL + 20F 2 + 4.632456F 3 - 13.264912F^ = 
-2 Fq - 4F X + 4F 2 - 0.632456f 3 + 2.632456F4 « S 1 
F 0 - 6 f 2 + 2.3l6228F 3 + 1.367544Fj^ = S 2 






( 71 ) 
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and 


[j?0 + 21.6 (oq ' + p 0 ' + 70 ' ) ^0 “ [32 + 21.6(00’ + Pq ' + 27 o')J w i 
+ (8 + 21.67 0 ')w 2 = p’ 

-(8 + 10.8p^')vo + [24 + 21.6(0-1* + Pi 1 + 7i’)>i 
■ jl6 + 21.6 (oi' + 7i')Jw 2 » p' 

2w q - Jl6 + 10.8(02* + Pg'jjwi + j^20 + 21.6(02' + P 2 ‘ + 72’)jw2 = 


A 72) 




It Is to "be noted that the terms of the left-hand side of equation (71 
do not change if the assumed values of K and S are changed. 

Equation (71) can he solved uniquely, therefore, in terms of K's 
and S's. The given, auxiliary, and final matrices obtained by Crout's 
method are given in tables 1, 2, and 3, respectively. More significant 
figures than required are used to ensure good results. 


The solutions of equation (71) are as follows: 


F 0 = -O.O487O3K0 
F 1 = -0.111203Kb 

F 2 = -0.103085KQ 

F3 = -0. 18993 7Ko 
F^ = -0.094968KQ 


O.265696K1 

0.307363Ki 

O.3H962K1 

O.506498K1 

0. 3165611^ 


O.225IHK2 

0.2355271^ 

0.2210522^ 

0.316561K2 

0.269077^ 


o.304n4Si 

0.262447S 1 

0.16288OS1 

0.253249S 1 

0.063312S1 


For a numerical example of the computation, let 


0.309525S 2 


0.288692S 2 

0.317642S 2 

0.126624S 2 

0.221593S 2 

J 


H73) 


p = 100 


p* = 0.0421875P = 4.218750 
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From the curve a for wq ~ 
it ie estimated that 




(figs. 11 and 12), 


w Q = 1.135 
w 

- - 0.7535 
vq 

w 2 

— = 0.5775 
w 0 

The first trial values are 

v Q = 1.135 

v ± = 0.855222 

v = 0.655463 
2 

These values are written at the right-hand corners ‘below the corre- 
sponding net points. The finite -difference patterns are used as given 
in figure 4, and a, 0, 7 , v n+1 - v n , and then K and S are 

found at the net points (fig. 13) • As an example, 

= -2(1.135000 - 0.855222) - -0.559556 
7 0 = 1.135000 + 0.655463 - 2 X 0.855222 = 0.080019 
Kq = (0.080019) 2 - (-0*559556) 2 = -0.306700 

Similarly, it is found that 

K-l = -0.189997 

Kg = 0.221966 

s. = 2.368276 
S 2 = 1.373368 



62 


MCA TN No. 1425 


From equation (73) the values 



of F ' a are obtained as follows : 

-1.129866 

-0.977802 

-0.780162 

-0.689444 

-0.424723 


These values are substituted in any one of the expressions (equation (71)) 
as a check and then are recorded at the net points, as in figure 13. 
Similarly, the values of a', P ' , and 7 ' are recorded below the 
corresponding values of F. 

Equation (72) can now be written and the given matrix is 




w. 


p » Check 

r column 


34.122771 

-12.269024 

2.000000 


-47.107213 

36.930948 

-19.408458 


8.984442 

•20.392900 

28.313451 


4.218750 

4.218750 

4.218750 


0.218750 

8.487774 

15.123743 


The check column can be obtained by using the following relation: 


Check column 


-4 + p* 

10.8p 1 * + p' 


6 + 21.6(02' + Pg') + p* 


The sian of the elements in a row should be equal to the value of the 
element of the same row in the check column. This procedure provides a 
check for the substitution made in the given matrix. 
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The first approximation gives, therefore 

w 0 = 1.117078 

= 0.843225 
w 2 = 0.648112 

A computation similar to the one outlined in the foregoing numerical 
example gives 


Kq = - 0.293781 
= -0.184115 
= 0.214841 
s 1 = 2.299072 

S 2 = 1.339974 


As a second trial, assume 

Kq = ^(-0.306700 - 0.293781) = -0.300241 
= |< -0.189997 - 0.184115) = -0.187056 

= i{0. 221966 + 0.214841) = 0.2181.04 

5 1 = 2 *368276 + 2.299072) = 2.333673 

5 2 = i(l. 373368 + 1.339974) = 1.356671 

The results of the second, third, and fourth trials are shovn 
in figure 13 . The corresponding assumed and computed values of the 
fourth trial are 
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Assumed 

Computed 

«0 

-0.300446 

-0.300006 

% 

-.187407 

-.187472 

K 2 

.218738 

.218786 

S 1 

2.3379^1 

2.338531 

S 2 

1.360090 

I.36063I 


The first three values check with one another, and the results, corrects* 
to the third decimal place, are 


v 0 = 1.1269 
w^ = 0.8502 
= O.6528 


The Large - deflections problem, n = 3 »~ When n is taksn to be 
greater than 2, the same procedure of computation as that in the case 
of n = 2 is still valid. As an example, the case of n = 3 will be 
considered, when the square plate is subjected to a uniform pressure 
of p = 100. 


After using the boundary conditions, the two sets of difference 
equations (6l) and (62) are obtained. Equation (6l) can be solved in 
terms of K’s and S's, and the results are given in table 4. 


From the curves of Wq ~ p. 


w 0 ~ P > 


wq 


P, 


^3 

w o 


v 4 

P> ~ 
w 0 


Vc 

and ~ p (figs. 12 and 14), the following values are obtained by 

"A 


extrapolation : 


Wq m 1.124? 


= 0.8891 

WQ 

Wp 

£ = 0.7932 
w 0 



0.5516 
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— = 0.5037 

w 0 


Y5 

w O 


0.3497 


For a first trial, it is assured that 


v 0 « 1.124700 
w-l - 0.999971 
« 0.892112 
w^ = O.620385 
V k - 0.566511 

v 5 - 0.393308 

Again these values are written at the right-hand corners "below the 
corresponding net points. With the computed values of a, 3, 7 , 

Ajv, and A^w, the following values are obtained: 


K<) = -0.061945 
» -0.052063 
Kg « -0.024186 
k 3 = -0.023043 

K4 = 0.001252 
= 0.106245 
S x = 1.592696 
s g = 1.282838 
s 3 = 0.548700 

By table 4 the values of F's are found to be 


F q = -1.095495 
F 1 = -I.O28996 
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Fg = 


-0.950911 
-0.868159 
- 0.762520 
-0.505761 
-O.67585O 
-0. 51+6620 
-0.239729 


The values of F's are written at the left-hand corners below the 
corresponding net points, and the values of a*, 0', and 7 1 are 

computed . 

When the values of a', 0', and 7 * are substituted into 
equation (62) and it is noted that p' = 0.00833333P = 0.833333, the 
given matrix of the equations is obtained as in table 5 and the auxiliary 
matrix as in table 6, and the solutions of equation (62) given by the 
final matrix are 


w 0 = 1.123384 
• 0.998956 
w 2 * 0.891465 
w^ = 0.620342 
v 4 = 0.565591 
v 5 « 0.390999 


It might be pointed out here that the check column of the given matrix 
may be obtained by a direct substitution by using the following relations: 
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Check column 
P* 

-1 + p* 

-2 + p' 

2 X 10 . 80^ ' + p* 

1 x 10.8P 5 ' + p' 

6 + 21.6(a 6 ' + 75') + p‘ 


This procedure vould provide a way of checking ‘the substitution in the 
given na.tr ix, since the sum of the elements in any row should he equal 
to the element of the same row in the check column. 

The values of Kq, K^, £ 3, Kjj., K5 , Sjl, S 2 , and are 

found from the computed values of w's. The mean values of the K's 
an d S's first assunBd and those computed are used as the trial values 
for the second cycle, and so on. At the end of the third trial, the 
following assumed and computed values are obtained: 



Assumed 

Computed 

*0 

-O.061763 

-0.06l695 

% 

-.051947 

-.051894 

*2 

-.024660 

-.024799 

*3 

-.023377 

-.023477 

K 4 

.001614 

.001697 

*5 

.106177 

.106204 

S 1 

1.592106 

1.592078 

S 2 

1.281878 

1.281814 


.546560 

.546173 
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These values check with one another to the fourth decimal place. The 
deflections at the various net points, accurate to the fourth dec imal 
place , are 


w q = 1.1240 
w x = 0.9995 
w 2 = 0.8920 
w^ = 0.6207 
W^ = 0.5660 
v 5 = 0.3915 


The results of various trials are shown in figure 15 • 


RELAXATION METHOD 


When a more accurate result is needed, the plate must he divided 
into a set of finer nets. The number of simultaneous equations 
increases as the number of nets is increased. In order to avoid the 
solution of simultaneous equations, Southwell's relaxation method 
may be used. The so-called relaxation method is essentially a clever 
scheme for guessing the solution of a system of difference equations. 

A brief description of the method and a numerical example, the small - 
deflection problem of a square plate, are given in appendix A. 

The solution of the general case of the large -deflection problems 
of rectangular plates by the relaxation method has been studied by Green 
and Southwell and their method was outlined previously. Green and 
Southwell worked with the three complicated equilibrium equations in 
terms of the displacements u and v and the deflection w. However, 
it is satisfactory to use the two much simpler equations in terms of the 
stress function F and the deflection w. 

The fundamental differential equations (l) and (2) can be rewritten 
as follows : 


vV = k (74) 

V 4 * = 10. 8p + 10.8k' 


(75) 
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where 


k = 


k' 


( 

b^v b^w 

\^x by) 

dx 2 dy 2 

d 2 F d£w 

d 2 F b^v 

dx 2 by^ 

C\J 

A 

OJ 

& 


d 2 F d 2 w 
bx by bx by 


In applying the relaxation method, as usual, the domain of the problem 
to he solved is first drawn, and the net points chosen. Since there 
are two simultaneous equations to he solved, two sheets of paper may he 
used, one for F and one for w. A set of solutions of F and w are 
guessed enH are recorded on the F~ and w “plane s , respectively. 

By starting from the assumed values of w, K can he computed 
without difficulty. Equation (74) is then a linear differential equation 
for F, and the hihamonic relaxation pattern may he used. After the 
residues at each point have heen reduced to the desired extent, the 
new values of F may he substituted into equation (75) end it may he 
solved by the relaxation method. Equation (75) leads to a rather 
complicated relaxation pattern for w. In actual computations the 
biharmonic pattern may he used, the assumed values of w being used for 
the computation of k'. By means of the relaxation process, the residues 
at all points are reduced somewhat. New values of k' are computed 
and the residues are then corrected. The relaxation operation is applied 
again until the values of w are determined to the desired accuracy. 

The average values of the new K's and S's and the originally assumed 
ones are now Used in the second cycle. The cycles are now repeated until 
the final results have the desired accuracy. 

In general, the boundary conditions for F are usually difficult 
to handl e. It is possible, however, to solve the boundary values of F 
in terms of its values for interior points. The boundary values of F 
vary from time to time as the interior valves change. The operation 
is rather complicated, but it can be handled. 

In the case of a square plate with given edge displacements, the 
bo undar y conditions as given by equation (49) may be used to some 
advantage. Equation (74) can be written as 

= k ] 


V 2 f = T 


( 76 ) 
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and the boundary conditions are given by 


and 


T 0 1 + 2T, . + . . . + 2T + T . = S. 

l,i m-l,i m,i i 


m 


(1 - p)(AZ) 2 


n=0 


w n,i - w n-l,i 


> 


(77) 



In using this form, not only are the boundary conditions much easier to 
handle, but also the relaxation pattern is simplified from the biharmonic 
type to the harmonic type. The simplification is obtained at the 
expense of introducing one more equation into the system and therefore 
considering one more plane. The results obtained are given in figure l6. 


DISCUSSION OF RESULTS 


The bending problem of a square plate under uniform normal pressure, 
with the edges prevented from displacements along the supports but free 
to rotate about them, is studied by the finite -difference approximations. 
The difference equations are solved by the method of successive approxi- 
mation and by the relaxation method. The computation starts with n = 1 
to n = 3 j in which case the plate is divided into 36 square n^ts with 

25 inner points. The maximum normal pressure calculated is - = 250. 

Eh 4 

After the values of w and F have 
can be found by the following relations: 

„ , & yfr 

1 3 y 2 (Al ) 2 

0 , _ V* 

y dx 2 (Al ) 2 


been determined, the stresses 


P* 


(Al)‘ 


(Al)' 
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a x “ " 


2(1 


+ 


2(1 - ^ 2 )(A l) 2 


(a + nP) 


y 2(1 - n 2 ) ( Al ) 2 


(P + pa) 


where a' and a" are the membrane stress and the extreme -fiber 
bending stress, respectively. The total stresses a are the sum of 
the nembrane and bending stresses at the section and are maxim u m at 
the extreme fiber of the plate . They are 


a x - °x‘ + a x" 


CXy = by' 4- by" 


At the center of the square plate, a* = P* and a = p, and therefore 
the stresses are 


, , a ' 0* 

° x m °* = (aT) 5 ’ (57* 


II _ If 


= 3 2(1 - 4)(A1) 2 = 2(1 - H)(A0 2 


The deflections at various points determined in the cases n = 1, 
n = 2 and n = 3 are tabulated in tables 7 to 9. The center deflections 
are plotted against the normal pressure ratio in figure 12. The membrane 
stresses in the center of the plate and at the centers of the edgss are 
tabulated in table 10 and are plotted in figure 17. The bending and 
total stresses are tabulated in table 11 and are plotted in figure 18 . 
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A study of the results shows that the maximum error In center 
deflections is 0.4-7 percent for n = 2 in comparison with n = 3 
and the maximum error in the center membrane stresses is 0.44 percent, 


"both values being conservative. Both maximum errors occurred at 

pa^ 


Eh^ 
= 12.5 


= 250 . 


The error in the center bending stresses is 2 percent at 

pa 4 Eh4 

and is O.83 percent at i —~ = 250, both values being unsafe. The. error 

Eh^ 4 

in the center extreme -fiber stresses is 1.6 percent at r = 12.5 

pa 4 Eh 

and 0.17 percent at • — = 250, both values being safe. The error in 

Eh 4 

the membrane stresses at the center of the sides is 12 percent for both 

4 

a 'a 2 /Eh 2 and a * a 2 /Sh 2 at = 12.5 and 8.9 percent for both 

7 Eh\ 

a x 'a 2 /Eh 2 and a 'a 2 /Eh 2 at = 250, these values being unsafe. 

* Eh 4 


4 

At ^ 
Eh* 


One case of n = 4 has been solved by the relaxation method. 


100, it is found that 


= 1 . 1250 , 


0 x ^* a2 
*0 

Eh 2 


J JQ 


•«2 


Eh* 


= 4.786, 


a 0 a 2 cj Xi 'a 2 a yi' a2 

= 11.394, 1 = 9*688, and — 1 = 3 .064. When the results 

Eh 2 Eh 2 Eh 2 4 

for n = 3 are compared with those for n = 4 at — r = 100, the 

Eh 4 


center deflection has an error of 0.09 percent, the center membrane stress 
has an error of 0.02 percent, the center total stress has an error of 
0.5 percent, and the membrane stresses a x 'a 2 /Eh.2 and o-y 'a 2 /Eh 2 have 
the errors of 4.2 percent and 4.1 percent, respectively, all values 
being unsafe. Since in the present case only the center deflections and 
stresses are to be investigated and the errors are sufficiently small 
for engineering purposes, the case n = 3 is considered to be satisfactory 
for the final results. 


The center deflections obtained by Way (reference 15). Levy 
(references 17 and 19) , and Head and Sechler (reference 23) are plotted 
in figure 19 for conqjarison with the present results . The center 
membrane , bending, and total stresses are plotted in figure 20 to 
compare with the results by Levy (references 17 and 19)* It is seen 
from these results that the center deflections are in good agreement with 
test results from the California Institute of Technology up to 


= 120 . 

Eh 4 


The theoretical results 3eem to be too low at higher pressures. 
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It is interesting to note that the test results are really for 
clamped -edge plates. The clamping effect seems to he only local, and 
at the center of the plate the plate behaves Just as though it were 
simply supported] that is, the plate is free to rotate about its edges. 

From the point of view of the engineer designing the plate, the 
total stresses at the center of the edges are still much larger in the 
case of clamped edges than in all the other cases] hence, a design based 
on those stresses would give a conservative structure. The center 
deflections, however, would give an idea of the magnitude of the 
washboarding of a boat bottom while a seaplane is taxying or landing. 


CONCLUSIONS 


The following conclusions may be drawn from a theoretical analysis 
of an initially flat, rectangular plate with large deflections under 
either normal pressure or combined normal pressure and side thrust: 

1. The large -deflection problems of rectangular plates can be 
solved approximately by the present method with any boundary conditions 
and to any degree of accuracy required . Although it is still difficult, 
the present n»thod is, nevertheless, simpler than the previously used 
methods for giving the same degree of accuracy. 

2. For the square plate considered, case n = 3 gives results of 
good accuracy, and the results are consistent with the existing theories. 

3. The clamping effect of a clamped thin plate seems to be only 
local. At the center, the plate behaves more like a plate with simply 
supported edges] that is, the thin plate is approximately free to rotate 
about its edges. 

pa.4 S pa 4 

4. The test results show that, at - ei r > 175 ( where — r is 

\Eh 4 V Eh 4 

nondimensional form for normal pressure J , ail. the existing solutions of 

the differential equations give unsafe results for center deflection 
for a square plate . This conclusion perhaps suggests the range in 
which the differential equations may be applied. 


5. The present results 
stresses give good agreement 


of the center deflections and membrane 

k 

P a 

with the test results when — r < 120. 

Eh 


Massachusetts Institute of Technology 
Cambridge, Mass., March 4, 1946 
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APPENDIX A 

BRIEF DESCRIPTION OF RELAXATION METHOD 


Tho idea behind the treatment by the relaxation method is 
essentially just the same as that by Cross' method of moment distribution 
in the case of bending of continuous beams . It seems, therefore, 
easiest to explain the relaxation method by a comparison with the 
moire nt -distribution method, since the latter is well accepted and is 
familiar to most structural engineers. 

The redundant beam as shown in figure 2l(a) is now examined. The 
procedure for obtaining the redundant support moments by the moment - 
distribution method is well known. The first step in the moment- 
distribution analysis is to assure that the slope at each of the four 
supports is zero. By this assumption, the end moments at A, B, 

C, and D can be found without difficulty. Tho result is shown 
in figure 21(b). Here the boundary conditions at A and B are 
satisfied, and the principle of continuity Is also satisfied. The 
condition of equilibrium, however, is not satisfied, since there are 
unbalanced moments at B and C. The moment -distribution method now 
offers a procedure to balance these unbalanced moments by a relaxation 
based on consistent deformations. The analysis by the relaxation 
irethod, in this case, would be essentially the same. The moments at A, 
B, C, and D are ass ume d to satisfy the boundary conditions and the 
condition of continuity. The unbalanced moments at B and C are 
then distributed by the relaxation based on consistent deformations. 

The difference lies in that the relaxation method offers more freedom 
in assuming the end moirents and therefore could make the convergence of 
the operations more rapid. On the other hand, however, it becomes 
difficult to assume these values. 

The method of moment distributions applies only to redundant 
structures, but the application of the relaxation method extends much 
further, and its application to the partial differential equations 
has brought the study of engineering sciences into a new era because 
the boundary conditions are now no longer difficult to be described and 
to be satisfied. 

The procedure can be Illustrated by a study of the small -deflection 

w * 

theory of thin plates. Letting w = where w' and p are the 

nondimensional deflection and pressure, respectively, gives the 
following equilibrium equation in terms of the finite difference 


A^w + 2A xy 2 w + = 12(1 - n 2 )(Al) 4 


(Al) 
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In order to solve the problem, the domain to be Investigated is drawn 
and the net points chosen. Values of w are assumed to satisfy the 
boundary conditions and are then written adjacent to each point of the 
net. From these values of w, the residuals Q at points (m,n) 
are computed and recorded as follows : 


Sn,n = 2( ^ w m,n " ®(^m+l,n + w m-l,n + w m,n+l + w m,n-l) 

+ 2 (^m+1 , n+1 + w m+l,n-l + W m-l,n+l + W m-l,n-l) 

+ (*m+2,n + v m-2,n + %,n+2 + w m,n-2) 

- 12(1 - n 2 )(Al) 4 (A2) 


The residuals Q thus computed can be thought of as an unbalanced 
force which must b« removed from the system. Now, instead of setting 
up a specific iteration process, it is merely observed that if the 
deflection at one point (m,n) is altered, all others remaining fixed, 
the residuals will change according to the pattern of figure 4, the 
relaxation pattern. Each change of w at any point effects a redistri- 
bution of the residuals Q among the net points, and such changes of w 
are de sired as will move al 1 the unbalanced force 0 to the houndary . 

For a simply supported plate, the deflection and bending moments 
are zero along the edges. Equation (Al) can be written as 

V 2 ^ 2 ^ = p 


Letting V^w = M makes possible the formulation of the boundary -value 
problem as follows: 



> 


M = 0 along the four edges J 


(A3) 


and 


'A 


^w = M 

w = 0 along the fora* edges 


(A4) 
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The problems can now be solved In two steps, that is, first, by 
use of equation (A3) and then by use of equation (A4) . This transformatio 
greatly reduces the labor required in applying the relaxation method 
because the relaxation pattern of the harmonic or Laplac ian type is 
much simpler than that of the biharmonic type. 

As an example , the boundary -value problem is solved when the plate 
is a square one. The process is considered with n = 4. From the 
previous results as found from the calculations with n = 3, the values 
of w at all the net points can be assumed. By equation (A4) 


^Vjn w m+l,n + w m-l,n + w m,n+l + w m,n-l " ^ w m,n 

The values of Mjjj n are then recorded at the right of the corresponding 
net point, and the residuals 

«m,n " *Vl,n + »m-l,n + »m,n + l + * ^^l) k (A6) 

are computed and are recorded at the left of these net points. The 
results are shown in figures 22(a) and 22(b). For example, 

Mq = 4w-^ - 4 wq = 4(0.04o6) - 4(0.0437) 

= -0.0124 


m 4 = w 2 + w 3 + w 5 + w 7 “ ^ w 4 


= 0.0377 + 0.0316 + 0.0231 + O.OI63 - 4(0.0295) 


= -0.0093 


Qq = 4m x - 4Mq - 0.002637 


= 4( -0.0117) - 4( -0.0124) - 0.002637 


= O.OOOI63 
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Q4 = Mg + + My - 4M4 - O.OO2637 

= -0.0106 - 0.0093 - 0.0078 - 0.0064 - 4( -0.0093) - 0.002637 

= 0.001463 

where O.OO2637 = 12(1 - ^ 2 )(Al)^, since n 2 = 0.1 and Al = ^ « 

The largest counterbalanced M occurs in the vicinity of the 
greatest deviation of the assumed values from the correct solution; 
so changes are first made at this point. An examination of figure 22(h) 
shows that the greatest residual occurs at point 2* Since 


Q 2 = 2M-L + 2M4 - 4M2 - 0.002637 


a change of Mg would change Qg "by an amount equal to four times 
(-AMg). Mathematically, 

AQ2 = - 4 AMg 

where A denotes the amount of change. Adding -0.0004 to Mg while 
assuming all the other values of M to remain unchanged gives 
AQg = 0.0016, and Qg is now equal to -0. 000637* If a nomenclature 
similar to that in the method of moment distribution is used, this 
process can "be called balancing the unbalanced Q. A symbol (bl) is put 
at the side of the value to indicate the first balancing. Now it is 
observed that 


Q x = Mq + 2 Mg + M 3 - 4M x - 0.002637 


and 


Q4 = Mg + M3 + M5 + My - 4M4 - O.OO2637 

A change of Mg with all the other M's fixed would change Qy_ and Q4 
by the relations as follows: 
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^4 = £M 2 


Now, "by relaxing the nets, 

= 2 (-0.0004) = -0.0008 
AQ4 = -0.0004 


and 


Q x = 0.001263 - 0.0008 = 0.000463 
= 0.001463 - o.ooo4 = 0.001063 

These operations may he called carrying-over and he denoted hy (cl). 

The whole process consists of 20 balancing and carrying-over 
operations hy similar calculations. The detailed operations of the 
computations are shoun in figure 22(h). After the values of M's are 
computed, the residuals are computed as follows: 


%i,n' “ w m+l,n + w m-l,n + W m,n+1 + W m,n-1 “ lw m,n " ^n 

The values of w may he determined hy a similar series of calculations 
The detailed operations and computations are shown in figure 22 (a). 

The whole process consists of 11 balancing and carrying-over operations 
The center deflection ratio thus obtained is, for n = 0.316228, 


v Q = 0.043790p 


For p = 0.3, 


w n = 0.043790 X ^2ip 
U 0.9 


= 0 .044 3p 


which checks exactly with the exact analytical solution. 
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For thin plates with clamped edges, the boundary conditions are 


v = 0 



•v 

o 

II 

m 

along 

X 

11 

ro'fH 

o 

II 

along 

v 3 ii 
* 2 


The relaxation pattern of the biharmonic type must be vised in this case. 
Although the pattern is more complicated, the process is essentially the 
same . 


After the essential idea of the relaxation method is grasped, other 
problems may be solved by rather obvious steps. It may be noted that 
no question of convergence can occur in the general relaxation process 
since no specific instructions are given. If, after some steps, the 
residuals get worse, the intelligent computer mates changes in the 
opposite direction. These remarks, however, oversimplify the problem 
socewhat because of two facts: first, the computer may become confused 

as to whether the residuals are really better, and, secondly, there is 
always a question of whether a solution with zero residuals exists. 
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TABLE 6 .- SOLUTIONS OF EQUATION ( 62 ); AUXILIARY MATRIX 
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TABLE 7." CENTER DEFLECTIONS 
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TABLE 8.- DEFLECTIONS AT VARIOUS POINTS 
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TABLE 9.- DEFLECTIONS AT VARIOUS POINTS 


[n = 3] 


pa^/Eh^ 

w 0 /h 

w x /h 

w 2 /h 

V3/I1 

v k /h 

v 5 /h 

0 

0 

0 

0 

0 

0 

0 

12.5 

.4055 

.3564 

.3136 

.2139 

.1890 

.1159 

25 

.6083 

.5365 

.4738 

.3249 

.2892 

.1822 

50 

.8460 

.7494 

.6650 

.4592 

.4131 

.2711 

75 

1.0031 

.8905 

• 7930 

.5500 

.4986 

.3370 

100 

1.1240 

• 9995 

.8920 

.6207 

.5660 

.3915 

150 

1.3104 

1.1677 

1.0450 

.7305 

.6717 

.4804 

200 

1.4557 

1.2988 

1.1641 

.8164 

.7551 

• 5531 

250 

1.5770 

1.4081 

1.2634 

.8880 

.8249 

.6149 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 




4 

Eh 


0 

12.5 

25 

50 

75 

100 

150 

200 

250 


a 't? 
% a 

% a 

a i a 2 

ff yi a 

a X!‘a 2 

Eh 2 

Eh 2 

Eh 2 

Eh 2 

n =» 2 

n = 3 

n = 2 

Mi 

KfeCl 

n » 2 

ra 

0 

0 

0 

0 

0 

0 

.6103 

.6089 

.3338 

•3795 

1.055 

1.200 

1.384 

1-377 

.7612 

.8574 

2.407 

2.711 

2.695 

2.683 

1.484 

1.66l 

4.693 

5.254 

3.806 

3.792 

2.096 

2.341 

6.628 

7.401 

4.802 

4.785 

2.643 

2.943 

8.357 

9.305 

6.566 

6.542 

3.613 

4.001 

11.43 

12.65 

8.136 

8.103 

1 

4.473 

4.929 

14.15 

15-59 

9.575 

9.533 

5.264 

5.778 

16.64 

18.27 




















92 


NACA TN No. 14-25 


TABLE 11.- EXTREME -FIBER BENDING AND 
TOTAL STRESSES AT CENTER OF PLATE 
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Figure 1.- Coordinate system and typical loading. Figure 2.- Coordinate system used by Levy 

and Greenman. 
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Figure 12.- Center deflections for a square plate under normal pressure 
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Figure 15.- Results of various approximations, p - 100; n = 3. 
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Figure 16.- Relaxation method. 
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Figure 17.- Membrane stresses. 
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Figure 19.- Center deflections as obtained by experimental and various theories. 
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Figure 20.- Center stresses as computed from various theories. 
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(a) Redundant beam. 
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(b) End moments for zero slope at each support. 
Figure 21.- Moment -distribution method. 
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Figure 22.- Solution of small -deflection theory by relaxation method. 



NACA TN No. 1425 


113 





NATIONAL 

. ADVISORY 


COMMITTEE FOR AERONAUTICS 

9 


. 0003*9 
(ho) . 000/43 
Orr) 400*03 

(C/j) 40000 / 

- 003 * 

. 0 OO*(ijf) 


-.009 000 

s 


* 


-. 00**97 
(bi)- ooo *97 

(c*)-jooc rrr 
(c*y-.oo<jorr 

(bef J 0 OO 0 O 7 
(C/ 4 ) 0 COOO/ 

Ccap-.ooo oom 



.000 7*3 

(<4) 

40 / 33 % 

(H) 

.000 /** 

Oa) 

.000 033 

(bn) 

.000 007 

(C/ 4 ) 

■0OOOO/ 

(a 9) 

3 

000003 



-4/0* .00/ 0*3 -.0099 -. ooo 797 -.0090 

-.OOOf 1%4> fc) .00/0*9 .000*3 (M &i)-.000 307 -.000/3 (34) 

-.ooo 00* (3 *rc*J .0003/3 -Q0 O 003 (3h) (34) .Ooo 033 -.000 000(1*) 

-.a// 004 <Zl: 0OO,eT -.404073 ** ""O'* ~: aotS3 . — 

fc*) .000/93 Ccsa} -.000 0/7 

(**) <000009 (%/*)-. 00000/ 

(C//J-.0O0 0O1 (C/j) . 000 oom 

(c/*)~ 4 O 0 0 O* (c *)-400 OO/ 

(CJ 1) -.000000 ( C /j ) -.000009 

(ct§)-4ooo/o (c*o)-.ooooo* 
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&3) .000033 400/3 (3») (ci) <000049 .00003/%*) (U)-. 000397 -.000*0/2*) (c*) .000*7* -i00400*(h») 
!?* *00 003 - OOP oo*(3 q)(U) .000 3/3 -.000003(3/7 ) (<*)-. 000077 -.ooooctT**) (CO) -.000007 . 000003(2/9) 
CcnJ-.oooooo C*a)-.ooo t07 - — — — (ce) .00**43 — (*+ .0000*3 ooooot/toj 

r*t)-.ooooo' fir# .oeov (t4-.oeooo ■,-■—•*** 'jr c/ „ .00004* JIT 

(<//) .00000* (3*) .00000/ •(c/4) <0004/3 

(c/*)-. 000 004 (cq) .000000 (3/3) 400004 

(3a) .000 004 (ex) -400 003 (kit) 400003 

(t /4) .000000 (cat- 330 003 

(too) -400 000 


(b) V 2 M = 12(1 - H 2 )( A ;) 4 s 0.002637; Q = v 2 M - 0.002637. 


Figure 22.- Concluded. 
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